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We present algorithms for the group independent reduction of group theory factors of Feyn- 
man diagrams. We also give formulas and values for a large number of group invariants in 
which the group theory factors are expressed. This includes formulas for various contractions 
of symmetric invariant tensors, formulas and algorithms for the computation of characters and 
generalized Dynkin indices and trace identities. Tables of all Dynkin indices for all exceptional 
algebras are presented, as well as all trace identities to order equal to the dual Coxeter number. 
Further results are available through efficient computer algorithms. 



CO 
CN 

o 

oo 

a: 



X 



i 



1 Introduction 



As the number of loops to which perturbative field theories are evaluated increases, the groupQ 
structure of the individual diagrams becomes more and more complicated. This problem has been 
recognized many years ago and on a group-by-group basis some very compact algorithms were 
proposed for their computation. Especially for the defining and the adjoint representations of 
the classical groups SU(N), SO(N) and Sp(N) these algorithms can be implemented rather easily 
in a symbolic program that will then give the color trace of a diagram as a function of the parameter 
N H . The disadvantage of these algorithms is however that these results give no information about 
group invariants and hence it is only possible for very simple diagrams to generalize the results such 
that they are valid for arbitrary groups and arbitrary representations. Hence a different type of 
algorithm is needed, if one would like a more general answer. That such information is useful can 
be seen from some recent calculations in QCD [|| in which the representation in terms of invariants 
could show immediately why extrapolations of lower orders in perturbation theory could not be 
successful. In addition the presentation in terms of group invariants is more general and needs 
hardly any new work when one needs to apply it for different groups or representations. The need 
for this kind of generality is clear, for example, from grand unification and string theory, where all 
semi-simple Lie groups may occur. 

We consider non-abelian gauge theories based on simple compact Lie groups. The extension 
to semi-simple algebras and additional U(l) factors is then straightforward. The gauge bosons 
are assumed to couple to matter in some irreducible representation R of the gauge group. The 
generalization to reducible representations is also straightforward. The group-theoretical quantities 
that appear in the initial expressions are the structure constants f abc (appearing in gauge self- 
couplings and ghost couplings) and the Lie-algebra generators T§ in the representation R, appearing 
in the coupling of the gauge bosons to matter. In this paper we consider only "vacuum bubbles" , i.e. 
diagrams without external lines. As far as the group theoretical factor is concerned, our results are 
relevant for any diagram whose external lines carry no gauge quantum number, or for the absolute 
value squared of any amplitude if one sums over the gauge quantum numbers of all external lines. 
The group theoretical factor of other diagrams can be obtained by multiplying the diagram by 
projection operators. 

The group theory factor of a vacuum bubble diagram consists of traces of a certain number 
of matrices TJ-;, whose indices are fully contracted among each other and with some combina- 
tion of structure constants. Our goal is to obtain an expression for this factor that is minimally 
representation- or group-dependent. In principle, this goal can be achieved as follows. 

1. Express the traces in terms of symmetrized traces. This can always be done at the expense 
of some additional factors f abc 

Now one may simplify the resulting expression further by observing that the structure con- 
stants can be viewed as representation matrices in the adjoint representation. This allows us 
to 

2. Eliminate all closed loops of structure constants f abc . 

This amounts to performing step 1. on traces of adjoint matrices TJ|. Step 2. can also be 
performed in an algorithmic way to arbitrary order. However, the algorithm is not identical 
to that of step 1 because of the special properties of the adjoint representation. 

1 Since we are dealing with perturbation theory we only encounter Lie algebras, and we are insensitive to the global 
properties of the Lie group. Nevertheless, following standard practice, we will often use the word "group" rather 
than "algebra". 
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3. Express the symmetrized traces in terms of a standard basis of symmetric invariant tensors. 
A Lie algebra of rank r has precisely r such tensors ||, ||] . 

At this point we have succeeded in expressing every group theory factor in terms of r + 1 
representation-independent quantities, namely the symmetric tensors and f abc . The repre- 
sentation dependence is encapsulated in terms of (generalized) indices ||. We show how 
these indices can be computed for any representation of any Lie-algebra to any desired or- 
der. This algorithm requires a convenient choice for the basis of tensors, which is not the 
mathematically more elegant "orthogonal" basis advocated in §. The result is also to a 
reasonable extent group-independent. The only way group-dependence enters is trough the 
(non)-existence of certain invariant tensors, but one may simply take all possible tensors into 
account, and only eliminate them at the end. The only problem is that the group SO {AN) 
has two distinct tensors of rank 2N. This case can rather easily be dealt with explicitly. 

Although our main goal has now been achieved, the result is expressed in terms of many 
combinations of symmetric tensors and structure constants that are not all independent. Un- 
fortunately there do not seem to exist many mathematical results regarding these invariants. 
In particular, we are not aware of any theorem regarding the minimal number of invariant 
combinations. For this reason the rest of our program is limited to finite orders, and is not 
guaranteed to yield the optimal answer in all cases. As a first step we 

4. Eliminate as much as possible the remaining structure constants f abc . 

We do not know of a proof that this is always possible, and in fact we have only been able 
to do this explicitly up to a certain order. The first object where we were unable to perform 
step 4 is built out of two structure constants and three rank 4 symmetric tensors. 

5. If step 4 is completed, one is left with a fully contracted combination of symmetric invariant 
tensors. We derive formulas expressing many such contractions in terms of a few basic ones. 

The last step is essentially group dependent (and therefore somewhat outside our main inter- 
est): 

6. Compute a formula for the basic invariants for each group. 

Here "basic invariant" is any of the contracted combinations of symmetric tensors that could 
not be expressed in terms of others, and any combination involving additional structure 
constants that could not be eliminated. 

None of these steps is new in itself, but we believe that in all cases we are going considerably 
beyond previous results (see e.g [7-25]). Since the application we have in mind is to Feynman 
diagrams, it is essential not just to develop an algorithm, but also to make sure it can be carried 
out efficiently. Complicated group theory factors appear only at higher orders in perturbation 
theory, which implies that one must be able to deal with a very large number of diagrams. 

The organization of this paper is as follows. In the next section we give some definitions and 
conventions, and present some well-known general results on invariant tensors. In section 3 we 
present the algorithm to perform step 1. Although this is in principle straightforward, without 
proper care such an algorithm may quickly get out of control. The same is true for step 2, which 
is presented in section 4. In section 5 we present the character method for computing indices and 
symmetrized traces. This section is based on results presented in |2(| and § , the main novelty 
being the extension to all higher indices of exceptional algebras. In section 6 we discuss steps 4 and 
5. Section 7 contains some remarks regarding the advantages and disadvantages of the two basis 
choices for the symmetric tensors. 
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In appendix A we present explicit results for indices and trace identities of the exceptional 
algebras; appendix B contains a description of the computer methods used, and in appendix C we 
give a few examples to demonstrate the efficiency of the algorithm. Appendix D contains many 
explicit formulas for invariants (step 6). In appendix E we discuss chiral traces in SO(2N) 

2 Generalities 
2.1 Definitions 

We consider simple Lie-algebras whose generators satisfy the commutation relation 

[T a ,T b ] =if abc T c (1) 

Our conventions is to use hermitean generators T a and to choose the Killing form proportional to 
d ab. 

Tr T a T b oc 5 ab (2) 

with a positive and representation dependent proportionality constant that will be fixed later. With 
this convention the structure constants f a b c are real and completely anti-symmetric. 
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Figure 1: Dynkin diagrams and labelling conventions 



Representations are denoted either by their dimensions (if no confusion is possible), or by their 
Dynkin labels (a\, . . . ,a r ), where r is the rank. Our labelling convention is indicated in fig |]. A 
bar denotes the complex conjugate representation; "A" denotes the adjoint representation, "i?" a 
generic representation and "r" the reference representation to be defined later. The dimensions of 
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these representations are denoted by Na, Nr or N r ; the quadratic Casimir eigenvalues (defined 
more precisely below) by Ca, Cr and C r respectively. 

The generators of the adjoint representation are related to the structure constants 



(T A ) a bc = ~i f abc • (3) 



2.2 Invariant Tensors 

We will encounter traces 

Tr ... T^" (4) 

in any representation R of any simple Lie algebra. We wish to express the result in the minimal 
number of quantities. 

Every trace defines an invariant tensor M: 

Tr T^ 1 . . . Tft n = M^- an . (5) 

This tensor is invariant because the trace is invariant under the replacement 

Tr - UrT^ 1 = UfT b R , (6) 

where Ur is an element of the group in the representation R; Ua is the same group element in the 
adjoint representation. Hence we have 

M ai...a„ = u a 1 b 1 _ _ _ jja n b n M b^...b n ^ ^ 
which implies, in infinitesimal form 

J2f cbai M^ 1 - h - an = , (8) 

i 

with b inserted at position i. This "generalized Jacobi-identity" may be taken as the definition of 
an invariant tensor. 



2.3 Casimir Operators 

Every invariant tensor M defines a Casimir operator C(M) 

C R {M)= T R ■■■T% n M ai - an . (9) 

ai,...,a n 

It follows from equation (||) that Cr(M) commutes with all generators Tr in the representation 
R. If R is irreducible Schur's lemma implies that Cr{M) is constant on the representation space 
R. Note that this is true independent of the symmetry of M, and irrespective of any concrete 
definition of M in terms of traces. All that is used is the Jacobi-identity (equation (j||)). 
Of special interest are the quadratic Casimir operators Cr, which we define as 

(T£T^ = CrSh . (10) 

As a special case of this identity we can write, using equation (|3|) 

jacdjbcd = CA6 ab _ 
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2.4 Symmetrized traces 

Not all invariant tensors and Casimir operators constructed so far are independent. We would like 
to express all traces in terms of a minimal set of invariant tensors. As a first step one may use the 
commutation relations to express the trace in a completely symmetric trace plus terms of lower 
order in the generators, which in their turn can also be expressed in terms of symmetrized traces. 
An efficient algorithm for doing this will be discussed in the next section. After this step we only 
need to consider symmetrized traces 

Str T ai T an = — V Tr T^ 1 ) T«M™) (12) 

7T 

where the sum is over all permutations of the indices (the cyclic permutation may of course be 
factored out using the cyclic property of the trace). 

For each representation one may define a symmetric invariant tensor with 

but this still vastly overparametrizes the problem, because a new tensor is defined for every order 
n and for every representation R. 



2.5 Basic Casimir invariants 

It is well-known that the number of independent symmetric invariant tensors is equal to the rank 
of the algebra. This can be seen as follows. For each invariant symmetric tensor d of order n define 
a polynomial 

P d {F) = F ai ... F an d ai - an , (14) 

where F a is a real vector of dimension equal to the dimension of the algebra, Na- The tensor d 
can be derived from Pd{F) by differentiating with respect to F: 



n 



]OF^-dF^ Pd{F) - (15) 



Although a priori Pd(F) is a polynomial in Na variables, the fact that d is an invariant tensor 
implies that P^ depends in fact only on r variables, where r is the rank of the algebra. This is true 
because the polynomial is invariant under 

pa _^ jjacpc ( ( 16 ) 

and it is well-known that for every F a one can find a transformation Ua that rotates F a into the 
Cartan subalgebra. Hence Pd{F) depends only on as many parameters as the dimension of the 
Cartan subalgebra, i.e. r. Therefore it is not surprising that any such polynomial can be expressed 
in terms of r basic ones, although the precise details (e.g. the orders of the basic polynomials) don't 
follow from this simple argument. 

The orders of the basic polynomials for each group are known [|], |[, and are given in the following 
table (for future purposes this table also gives the "dual Coxeter number" g). As explained above, 
each basic polynomial corresponds to an invariant tensor, which in its turn corresponds to a Casimir 
invariant. The implication of table [l] is that for any given algebra the polynomials 

Tr F% = Tr F a T^) n (17) 
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10, 


12,14,18 


E$ 


30 


2,8,12,14,18,20,24,30 



Table 1: Ranks of basic invariant tensors 



can be expressed in terms of r basic polynomials of degrees as indicated above. In section 5 we will 
show how to obtain such expressions for any irreducible representation (irrep) of any (semi)-simple 
Lie algebra. The ranks of the invariant tensors - or more accurately those^] that can be written as 
traces over some representation R - are in fact an outcome of these calculations. 

If we can express a polynomial corresponding to some invariant tensor d in terms of basic 
polynomials, we can also express the invariant tensors into basic ones. Namely, suppose 

P d (F) = T,U P ^F) (18) 

i 

where the sum is over various terms of this type, with coefficients. Then the differentiation (eqjl~5|) 
of a term on the right hand side yields precisely the fully symmetrized combination of the tensors 
di, with weight 1; this means that the overall combinatorial factor equals the number of terms. For 
example 

J^(£(n 2 ) 2 = l^r* + ^ bd + ^ bc ) • ( 19 ) 

e 

In practice we will express all higher traces in terms of r basic ones, but we do not obtain 
the full dependence on F of the basic traces, and consequently we cannot say anything about the 
explicit form of the basic invariant tensors. Given an explicit basis for the Lie- algebra one may 
compute the full -F-dependence, but that is the same as computing the invariant tensor directly by 
computing a trace. 



2.6 Indices 

Since a Casimir operator is constant on a irrep, its value can be computed by taking the symmetrized 
trace over this irrep. We will see explicitly how to expand a trace in terms of fundamental symmetric 
invariant tensors. In general one has 

Str T% . . . Tft™ = I n (R)d ai " an + products of lower orders . (20) 

The invariant tensors can be chosen in some representation-independent way, for example by com- 
puting it for one given reference representation. Then all symmetrized traces can be expressed 
in terms of this basis of tensors. The leading term necessarily has the indicated form, with a 
computable representation dependent coefficient I n (R). This coefficient is called the n th index of 

2 We are not aware of others, but also not of a proof that they do not exist. 
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the representation. If there is no fundamental invariant tensor of order n the indices I n {R) are 
obviously zero for any representation^] 

The extra terms in equation ( |20| ) are symmetrized products of lower order tensors such that 
the total order is n, without contracted indices. The coefficients of these terms will be called 
sub-indices. 

There is a lot of freedom in defining d ai "' an since we could have modified it by any combination 
of lower order terms in equation (|20|). Note that modifying the tensors by lower order terms does 
not affect the indices, but does change the sub-indices. This freedom can be used to impose the 
conditions g || 

d a,...a l ...a nd cn...a l = Q l<n (21) 

This then defines the symmetrized tensors up to an overall normalization. The normalization can 
be fixed by fixing a normalization for the indices. This basis will be referred to as the orthogonal 
basis. It is the most elegant one from a mathematical point of view, but, as we will see, not the most 
convenient one for our purposes. In the following we will use the notation (f^-— a i— an f or tensors in 
the orthogonal basis. 

As mentioned before, tensors defined in any basis can be used to define Casimir invariants, but 
using the orthogonal basis has a clear advantage because it leads to a simple relation with the 
indices (Nr is the dimension of R) : 

I n (R) = ^C P (R) , (22) 

where 

An = d a £- an d a £- an . (23) 

This relation holds because when contracting with an orthogonal tensor only the leading terms 
survive. Note that this is true even if we do not expand equation (|2^) in terms of the orthogonal 
basis, but in terms of any other basis. Hence the Casimir eigenvalues are determined up to a 
representation independent factor once the indices are known. 

The indices are of interest in their own right, as was in particular emphasized in ||. In some 
cases they have a topological interpretation via index theorems. Furthermore they satisfy a useful 
tensor product sum rule. If 

R x <g) R x = eRi (24) 

i 

then 

N Rl I p (R 2 ) + Nr 2 I p (Rx) =Y^Ip(Ri) ■ (25) 

i 

In subgroup embeddings H C G there is also such a sum rule for branching rules: if R — > ® r « 
then 

I^{R)=I p (G/H)Y^^{n) , (26) 

i 

where I p (G/H) is the embedding index (here we assume that both G and H have precisely one 
index of order p; other cases require just slightly more discussion). 

3 For D r , r even, there are two indices of order r. The additional one will be denoted as 7 r . We will deal with this 
case in more detail below. 
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2.7 Normalization 



To arrive at a universal normalization we make use of the following general formula for the quadratic 
Casimir invariant 

c r = 1 ED^ 2 )^ ( 27 ) 
i=i j=i 

Here are the Dynkin labels of the representation R, and Gij is the inverse Cartan matrix. The 
advantage of this formula is that the Dynkin labels as well as the Cartan matrix have a fixed 
normalization that is not subject to conventions. Only the overall normalization is convention 
dependent. The factor r/ is introduced to allow the reader to fix the normalization according to 
taste. The dependence on r\ will be shown explicitly in all formulas. Given the universality of (p7|) 
it is natural to choose r\ in a group- independent way. Using (|27]) get for Ca'- 

Ca = V9 , 

where g is the dual Coxeter number. 

This convention defines the normalization of the generators once we have fixed the rank 2 
symmetric tensor. The natural definition is 

d ab = d ab = § ah_ 



Then Af 2 = dfdf = N A - Now the second index I 2 (R) is also fixed via (||): 

HR) = . 

For the vector representations V of the classical Lie algebras we find then I 2 (V) = \r\ for SU(N) 
and Sp(N), and I 2 (V) = V for SO(N). 

There are (at least) two considerations that might lead to a choice for rj. First of all it is possible 
to fix the conventions in such a way that I 2 (R) is always an integer. This leads to the choice rj = 2. 
On the other hand there are standard choices for the generators of SU{2) namely T a = ^a a (where 
a a are the Pauli-matrices) , and for SO(N), namely T^ u = i(5^5j — 5j5"), where the pair [iv with 
[i < v represents an adjoint index. Unfortunately these two choices correspond to different values 
of r], namely r/ = 1 for SU(2) and r\ = 2 for SO(N). 

2.8 Indices versus Casimir invariants 

We conclude this section with a few historical remarks. 

A vast amount of literature exists on the computation and properties of Casimir invariants. 
Most of these papers, [7-19], give more or less explicit expressions for the Casimir eigenvalues of 



the classical Lie Algebras A n ,B n , C n and D n and in one case, |17], also for G 2 - In formulas 



for G 2 and F4 are obtained, whereas S§ was considered, up to order 14, in [25]. The issue of 
completeness of a set of Casimir operators was studied in 2C , 21]. 

In applications to Feynman diagrams indices are more important than Casimir invariants, be- 
cause traces over matter loops yield indices and sub-indices, and not Casimir invariants. Indices 
have been discussed most frequently in relation to chiral anomalies. The second index was intro- 
duced by Dynkin, [27], and generalized to higher order in |28[| . Shortly afterwards Q it was realized 



that the definition of the indices could be improved by imposing the orthogonality constraint (^). 
In H formulas are given for the indices of classical Lie algebras. Indices of exceptional algebras have 
been studied up to sixth order, mainly for the purpose of anomaly cancellation in ten dimensions, 
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relevant for string theory. Sub-indices, when defined in the orthogonal basis, can be expressed in 
terms of indices. Unfortunately these relations are difficult to obtain, and become very complicated 
at higher orders unless some lower indices vanish. In || formulas for SU (n) have been given up to 
fifth order. We have computed the sixth order formula, but the result is rather awkward and does 
not encourage extension to higher orders. 

Although, as explained above, indices are closely related to Casimir invariants, the available 
formulas for the latter are of little use to us since they do not use orthogonal tensors for the 
definition of the Casimirs. Even if they did, one would still need the normalization factor J\f n , (p3|). 
The computation of this factor for all Lie-groups and all values of n is a difficult problem, related 
to the even more difficult problem of determining the tensors (f^--- an explicitly. For recent progress 
on the latter problem for the classical Lie algebras see We do not present explicit expressions 
for the symmetric tensors here. Since they always appear in contracted form, we never need them 
explicitly. 

Furthermore the Casimir eigenvalues give no information on sub-indices. 

3 Reduction to symmetrized traces 

In this section we will discuss the reduction of traces as they occur in Feynman diagrams into 
the invariants of the previous section. This is by no means a trivial affair because the necessary 
symmetrizations make that the algorithms typically involve 0{n\) terms when there are n generators 
in the trace. It is therefore important to choose the method carefully. We will have to distinguish 
two cases. In this section we will make the reduction of traces of the type Str T^ 1 . . . T^ n in 
which R can be any representation with the exception of the adjoint representation. In the next 
section we will consider such traces for the adjoint representation. The special role of the adjoint 
representation lies in the fact that, because of the equation (|3|), the commutation relation 

[T%,T b A ]=i f abc T% (28) 

does not really diminish the number of generators of the adjoint representation. It is actually 
just a different way of writing the Jacobi identity. A related reason for considering the adjoint 
representation separately is that the reduction of the other traces generates new structure constants. 

We will continuously keep in mind that the algorithms we derive are for implementation in 
a symbolic computer program. This means that in many cases a recursion type algorithm may 
suffice, even though it may not be very practical for hand calculations. 

3.1 First stage elimination 

The first part of the reduction is dedicated to the replacement of the traces over the generators Tr 
by the invariants cIr. For all representations except for the adjoint this can be done in a general 
algorithm. One should realize however that for very complicated traces the results may not be very 
short. 

In general a trace is not symmetrized. Therefore the introduction of the tensors cZr needs some 
work with commutation relations to make it symmetric. On the other hand, computer algebra needs 
algorithms that work from a formula, rather than towards one. Hence one can use the substitution 

Tr[Tg • • • Tg\ = Tr[Tg ■ ■ ■ T* n \ - Str T^ 1 • • • T^ n + d%'" an (29) 

Writing out the symmetrized trace will of course give n\ terms, each with a factor \jn\. Then we 
can commute the various T^T^ till they are all in the order of the original trace after which the n\ 
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terms with n generators will cancel the original trace. At this point we are left with the symmetric 
tensor da and 0(n n\) terms which all have n — 1 generators. As a recursion it will eventually 
result in terms with only two generators for which we know the trace. This algorithm is however 
rather costly when the number of generators inside the trace is large. 

The above formula has as its main benefit that it proves that one can express a trace of 
generators T of any representation R ^ A in terms of symmetrized traces and structure constants 
/. For practical purposes we have a better algorithm. It is based on the formula: 

n /-i \j 

T^'---T^ } T h R = J2^ L B j E* 1 - anb (30) 

3=0 J - 

in which Bj is the j'-th Bernoulli number. ( Bq = 1, B\ = —1/2, B2 = 1/6 and B4 = —1/30) and 
the function Ej is defined by the recursion 

j^ai—On rp{ai r^a n } (SI) 

Ft Ft \ J 

n 

j^ai-- a„b ^""^ j^ai—a-i-iai+i—ank ^jaibk (32) 

i=l 

Basically in the Ej one extracts a string of j structure constants /. By writing out the functions 
Ej one can show with some work that the proof of this formula is equivalent to proving the relation 



Ra 



1 



n + 



i=0 



+ E r ■ ^H) 3 -,; , h .M + i-VfRn-i) (33) 

^(n + l-*)! %\ p* 

This equation should hold for any positive value of n and any choice of the Ri (The Ri represent 
the symmetric combination of T n to T' n with a T b inserted at i places from the right. Hence Rq 
is the left hand side term of equation (|30|)). Actually the Ri are independent objects and hence we 
have a set of equations each of which is characterized by a value for n and the index j of Rj. To 
prove the whole formula we have to prove that all of these equations are valid. We use the following 
approach: One can test their validity for any small value of n and all allowed values of j (computer 
algebra lets one check this easily up to n = 100). Next one takes the special case of only Rq not 
equal to zero and the case of only R n not equal to zero. These are rather easy to prove. Then one 
takes the case for R\ not equal to zero. This is only a little bit more complicated. Finally one can 
express the case for other values of n and j in terms of the equation for n — l,j — 1 and n, j — 1. 
This then combines into a proper induction proof. The above formula can give some interesting 
summations when one selects special values for the Ri like Ri = x 1 after which the inner sum can 
be done. 

An important part in the application of equation ( |30|) is how to terminate the recursion. We 
note that due to the cyclic property of traces 



Tr 



np\0,± rj-\ fin} rj-\b 



dp}-'^ (34) 



Additionally we can terminate two cases in which there are still two generators outside the sym- 
metrization: 

TiiTpfTpfTpf] = d a p} a2a3 + % -f aia2a *I 2 {R) (35) 
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Tr 



rp{oi rp 0.2} rpU3 rpO,4 

R R R R 



^0,10,20,30,4 _|_ ja-ia^k ( ja\a2k 



2' 

+ J_/ 2 (^)(y a i a 3 fc ja 2 a 4 k + jaia 4 k ja 2 a 3 fc-j ^gg^ 

All together this algorithm is far superior. For a trace of 7 generators T R all with different indices 
it is about 50 times faster than the one based on equation (|29|). 

Just as with traces of 7-matrices one can make the algorithms much faster with a number of 
supporting tricks for contracted indices. This is however not quite as simple as with 7-matrices. 
When the contracted indices are close to each other, we use 

T b R T a R T a R T R = C R T b R T R (37) 
T a R T b R T a R = (C R - l -C A )T b R . (38) 

The relation between ^(-R) and Cr is given by NrCr = I2(R)Na- Additionally we have 

rj-tCLl rj-t(l2 rpb rpC rpa\ rj-iCL2 rp&l rpQft rpb rpC rp(12 rp^l —("} r J nC / 

R R R R R R R R R R R R ^ R R R R \ J 

And then for contracted indices that are not very close to each other we can use: 



T R ..-T a R T b R T R = T R - ■ ■ T R T R T R + if b ^ c T R ■ ■ ■ T R T R 

= T R ... T R T R T R + if* c T R ■ ■ ■ T c R T a R - f^f^T R ■ ■ ■ (40) 

We commute the T R matrix also towards the T R matrix because 



= -\f h3c r d T d R 

= -\c A T b R (41) 

Hence we can always eliminate two generators T R when there is a pair of contracted indices. Some 
terms may however obtain two /-matrices in exchange. The above equation shows also that the 
trace of a string of generators of which two generators are contracted with the same structure 
constant /, will lead to simplifications by commuting the two generators towards each other. In 
this case many of the terms that come from a commutator have just one generator fewer than the 
original term. 

By now it should be clear why the adjoint representation cannot be treated in exactly the same 
way. Eliminating two generators at the cost of introducing two structure constants / leaves us with 
exactly the same number of generators of the adjoint representation. 

We can use a few extra shortcuts for simple cases to avoid the use of equation ( |30|) in those 
cases: 

Hr[TgT%Tg] = 4 iaaa3 + l -h{R)f aia2a3 

r Yl\T ai T a2 T CL3 T a4 ] — J a i a 2i3 a 4 _j_ \_ ( J0i04n ra2azn _ i020zn raia,4n\ 
L R R R R J R 2 ^ R * ) 

+— I2(R) (f aia4n f a2asn — y a i a 2« £030411^ (42) 

6 
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To get the last equation into its minimal form we have used the Jacobi identities: 



q _ yu*2iy«3Mj _|_ jhisj jiliij _|_ J*3Uj J«2Mi 
Q _ ^n«2iy«3«4j _|_ ^«2«3jy*l«4j _|_ ^«3«lj f'hUj 

It is possible to create similar shortcuts for the higher traces. This serves however not much 
purpose. The majority of cases involves short traces and these expressions are rather lengthy. 



4 Reduction of adjoint traces 

At this stage we have only tensors of the type d,R and structure constants / left as objects with 
indices. All these indices are indices in the adjoint space and hence all have Na dimensions. 
Additionally there can be various constants like the second order Casimir's and the second order 
indices hiR), but these do not play a role in the following. 

For the other representations the loops that define the trace were rather easy to find. For the 
adjoint representation this is more complicated: all three indices of the structure constant / can 
play a role and hence there are more possibilities. The advantage is that very often one can find 
'smaller' loops. If for instance we have a diagram that consists of only vertices in the adjoint 
representation and there are no loose ends we have the results of table ^ in which 'girth' is the size 



girth 


2 


3 


4 


5 


6 


7 


8 


9 


^min 


2 


4 


6 


10 


14 


24 


30 


54 



Table 2: Minimum number of vertices needed for a diagram with a given smallest loop. 



of the smallest loop in the diagram and n m i n is the minimal number of vertices needed to construct 
such a diagram. Of course in mixed diagrams in which other representations are involved one can 
get loops in the adjoint representation that have up to n/2 vertices if n is the total number of vertices 
in the diagram. Such would be the case of there is one loop of n/2 vertices in a representation R 
and a parallel loop of adjoint vertices. But in that case the loop is easy to find, and the symmetry 
of the invariant cIr that is present already makes the introduction of the dA into a triviality: due 
to its contraction with dft the trace over the adjoint generators has already been symmetrized. 

Let us first have a look at the canonical reduction algorithm of equation (|29| ) . At first one might 
be worried that it will not terminate for the adjoint representation. After all it does not diminish 
the number of structure constants /. One can see quickly however that the commutator terms have 
a simpler loop structure. Hence each term will end up having a number of /'s grouped into an 
invariant or have a loop with fewer /'s even though the total number of /'s is still the same. Once 
we have a loop with at most three structure constants it can be reduced with the equation 

ji\iia\ ji2ha,2 ji-shaz _ }_(J^faia20-3 (44) 

which can be derived from the Jacobi identity. Hence in all cases the number of structure constants 
/ will become less. Similarly the algorithm of equation (^) will reduce a number of generators of 
the adjoint representation to a symmetrized trace and a smaller number of structure constants /. 
Hence also this algorithm will terminate. 

It is however possible to be more efficient about the reductions inside the adjoint representation. 
Our first observation is that for loops with an odd number of vertices in the adjoint representation 
reverting the order of the vertices gives a minus sign. Hence the fully symmetric object with an 
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odd number of indices must be zero. And because of this we do not need a full symmetrization to 
express loops with an odd number of indices. For example for five indices we write 



^7i«ll5*4«3*2 

_ _piii2iaiAi5 _|_ pki^i^ pii$k 

j^piliikis ji-iizk _|_ piiizkic, ji^k _|_ pixiiki*, ji^k (45) 

in which F represents the trace over a number of f's. Its difference with a trace over the Ta- 
generators (||) is just powers of i. Similarly one can use the reversal symmetry for traces of an even 
number of /'s to reduce the amount of work by a factor two. 

Because of the above simplifications the efficiency of the two algorithms ((|2^) and (|30|)) is not 
very different for the adjoint representation. 

The reduction algorithm should be clear now. One looks for the smallest loop among the 
various /'s. Such a loop will not have contracted indices, because otherwise there would be a 
smaller loop. Hence we do not have to worry about contracted indices as we had to do for the other 
representations. If the loop has only two or three _f's, we can eliminate it with either equation (pd|) 
or equatio n (|44| ) . Otherwise we can use a simplified version of the canonical reduction algorithm of 
equation (^9|) to obtain an invariant and terms with a simpler loop. Actually the fastest way here 
is to tabulate this reduction all the way up to loops with 7 /'s. For loops of 8 or more vertices in 
the adjoint representation we use an adapted version of equation (|30|). 



5 Computation of symmetrized traces 

At this point our group theory factors consist of combinations of structure constants, symmetrized 
traces d/''"'° n over the adjoint representation, and symmetrized traces over one or more other 
irreducible representations. We will now show how such traces can be expressed in terms of r 
traces over a single representation, where r is the rank of the algebra. 

As explained in section 2, in principle there are three quantities one might be interested in: 
Casimir invariants, indices and symmetric tensors. 

The results presented here amount to a computation of the coefficients of combinations of 
fundamental traces appearing in the expansion of a trace in an arbitrary representation. In other 
words, we compute indices and sub-indices (but, as explained in section 2, the latter are basis-choice 
dependent). 

With our method these quantities can rather straightforwardly be computed to any desired 
order, and for any representation of any Lie algebra. To demonstrate this we will compute all the 
indices for the lowest-dimensional representations of the exceptional algebras, including the 30 th 
index of E%. 

The method we follow here is an extension of results of |2(| (where it was used to obtain the 
"elliptic genus" in string theory), which in its turn was an extension of results presented in ||| 
(where it was used for computing the indices of the classical algebras). 



5.1 Characters 

An extremely useful tool for computing traces are the characters 

Ch K (F) = Tr e FR (46) 
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where F R = F a T^. Hence the expansion of the exponential gives us all symmetrized traces in terms 
of the polynomials defined in equation (14). What makes the characters especially useful is their 
tensor property 

Ch Rl ® Rl (F) = Ch Rl (F)Ch R2 (F) , (47) 

which follows directly from its definition. In addition characters are combinations of traces and 
therefore also have nice properties on direct sums 

Ch Rl(SRl (F) = Ch Rl (F) + Ch R2 (F) . (48) 

With a little more effort one can also derive a formula for characters of symmetrized and anti- 
symmetrized tensor products |pS| |. These formulas can be derived from the following generating 
functions 

oo oo 

^x k Ch [khR (F) = det(l + xe F R ) = g exp(-(-x) / Ch i? (^)) (49) 

k=0 1=1 

oo oo 

x k Ch {khR (F) = det(l - xe F R )- 1 = JJ exp({x) l Ch R (lF)) . (50) 

k=0 1=1 

Here [k] denotes the order k anti-symmetric tensor product of some representation R, and (k) the 
order k symmetric product. We use the notation [k] * R or (k) * R to denote the anti-symmetrized 
or symmetrized tensor product of the representation R. Note that the sum in equation (|49| ) is in 
fact always finite. 

The generating functions can be expanded explicitly to obtain 



mil \ m 



CV)*fi= E il^T^ ^ J - (52) 



where the sum is over all partitions of the integer k into different integers rii, each appearing with 
multiplicity mj. 



5.2 Character computation method 

Our method for computing the characters is as follows. We begin by choosing a reference repre- 
sentation, which in all cases is the one of smallest dimension. The reference representations we 
choose for the simple Lie algebras are shown in table ^. (the last column of this table is explained 
later) Note that for SO(N) the reference representation is the vector representation for N > 7, 
but for lower values of N it is a spinor representation. Another way of saying this is that we treat 
SO(N),N < 6 according to the Lie-algebra isomorphisms D$ ~ ^3, B2 ~ C2 and B\ ~ A\. The 
last column of table § is discussed below. The algebras B%,Bi,Di and D5 are listed separately 
because, although they have the "standard" reference representation, they have non-standard index 
normalizations . 

For the reference representation the character is left in the form 

Ch(F) = Tr e F , (53) 

All traces whose order does not appear in table [l] can be expressed in terms of lower traces. Hence 
equation (|^) must be supplemented by trace identities for those traces. These trace identities will 
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Algebra 


Reference representation 


Dimension 


Indices 




(1,0,... ,0) 


r + 1 


1,...,1 


Bi 


(1) 


2 


1 


B 2 


(0,1) 


4 


1,1 


B 3 


(1,0,0) 


7 


2,2,1 


B 4 


(1,0,0,0) 


9 


2,1,1,2 


B r , r > 5 


(1,0,... ,0) 


2r + 1 


2,1,-.. ,1 


Cf 


(1,0,... ,0) 


2r 


1,-..,1 


D 3 


(0,0,1) 


4 


1,1,1 


D 4 


(1,0,0,0) 


8 


2,2,1,0 


D 5 


(1,0,0,0,0) 


10 


2,1,1,2,0 


D r , r > 3 


(1,0,... ,0) 


2r 


2,1,. ..,1,0 


G 2 


(0,1) 


7 


2,1 


Fa 


(0,0,0,1) 


26 


6,1,1,1 


Eq 


(1,0,0,0,0,0) 


27 


6,1,1,1,1,1 


E 7 


(0,0,0,0,0,1,0) 


56 


12,1,1,1,1,29,1229 


Es 


(1,0,0,0,0,0,0,0) 


248 


60,1,1,1,1,41,199,61 



Table 3: Reference representations, dimensions and indices. 



be derived below. The remaining traces will be called "fundamental" and equation (^) is taken to 
be the definition of the corresponding polynomials and symmetric tensors. This then defines a set 
of reference tensors: 

d ai...,a n = g tr T ax T a n ( 54 ) 

where n is the order of a fundamental Casimir operator. The precise form of this tensor, or 
equivalently the precise form of the fundamental polynomials Tr F n depends on the details of the 
Lie algebra basis choice, but will never be needed. 

Any other character is now written in terms of traces of F n over the reference representation, 
using all available trace identities. By differentiating with respect to F (c./. eq. ([0])) one can then 
read off the expression of any d,R in terms of reference tensors. 

This fails if the reference representation has an index that is zero. This happens only for the 
n th index of the algebra D n , and we will deal with that case separately. 

The next step is to express the characters of all "basic" representations in terms of the reference 
character. The i th basic representation is defined by Dynkin labels aj = 5{j, j = l,...,r. The 
most important tool for obtaining these characters is equation (|5l|). This yields all fundamental 
representations of the algebras A r and C r , whereas for the orthogonal groups only the spinor 
representations are still missing. The spinors, as well as the basic representations of the exceptional 
algebras, require some extra work, and are discussed below. 

Finally one can compute the characters of all other representations by using in a systematic 
way the sum rule for tensor products. It can be proved that for any simple Lie-algebra this allows 
one to relate the characters of all other irreps linearly to those of the basic ones. In principle this 
still allows for the possibility that complicated linear equations need to be solved. We find however, 
that one can organize the tensor products in such a way that only one unknown character appears 
at every step. This can be proved for the classical Lie algebras (see below), and we have checked 
it empirically for the exceptional ones. 
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Let us contrast this procedure with the computation via Weyl's character formula 

where A is the highest weight of a representation, the summation is over all elements w in the Weyl 
group W, e w is the determinant of w, p is the Weyl vector (with Dynkin labels all equal to 1), and 
h is a vector in weight space, which plays the role of F in the foregoing discussion. One obvious 
disadvantage of this formula is the summation over all elements of the Weyl group, although this 
is still manageable in most cases of interest. A less obvious disadvantage is that numerator and 
denominator both have a zero of order N + , the number of positive roots, in h. For example, to 
obtain the highest non-trivial Casimir eigenvalue of E%, which is of order 30, one needs to expand 
numerator and denominator to order N + + 30 = 150. This is an impossible task. The method 
sketched above, and worked out below, does allow an expansion of the character to order 30, even 
for Eg. 

An important ingredient in our procedure is obviously the computation of tensor products. 
Conceptually this is certainly not easier than the computation of characters, but nowadays computer 
programs exist that can do this very efficiently^ Rather than using characters (and in particular 
index sum rules) to compute tensor products, it is then more efficient to use tensor products 
to compute characters. The procedure described here requires just a small effort to compute 
the characters of the basic representations up to a certain desired order. The computation of 
the character of any other representation is then just a matter of simply polynomial operations 
(multiplications, additions and subtractions which can be efficiently performed by any symbolic 
manipulation program, such as FORM) guided by the output of a program that computes tensor 
products| 

We will now discuss the various types of algebras in more detail. 
5.3 A characters 

Let us now apply these tools first of all to Lie algebras of type A r (SU(r + 1)). For the reference 
representation we choose the vector representation (r + 1). Using equation ( |49| ) we can then 
immediately write down the characters for all the anti-symmetric tensor product representations 
[k]. In terms of Dynkin labels these are all the representations with labels (0, . . . , 0, 1, 0, . . . , 0), i.e. 
a single entry 1. These are precisely the basic representations. 



Now we can systematically use the tensor product rule (47) and the sum rule ( jig ) to obtain 
character formulas for all other irreducible representations. If 

Ri®R 2 = J2 eriiRi (56) 

i 

then 

Ch Rl (F)Ch R2 (F) = niCh^ (F) (57) 

i 

By computing the product of two known characters and subtracting the known characters on the 
right hand side one is left with the character of some (in general reducible) representation, which 
is thereby determined. 

4 We have used the programs LiE [^o[ that computes tensor products of Lie-algebra representations directly and 
Kac Q 

that uses the Verlinde formula to compute fusion rules of Kac-Moody algebras. A subset of these fusion rules 
coincides with tensor product rules, and it turns out that this precisely includes the tensor product rules we need. 

5 We have implemented this idea in the program Kac. The results in the appendix were produced in that way. At 
a given Kac-Moody level, the fusion rules that coincide with tensor products are found to be sufficient to o btain the 



characters of all representations at that level. For examples and software see http:/ /norma. nikhef,nl/~t5£ 
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To show how this works we label the SU(r + 1) irreps by Young tableaux and assign a partial 
ordering to them. We use Young tableaux because for A r they provide a convenient description of 
the tensor product rule. A Young tableau is ordered above another one if it has more columns; if 
the number of columns is the same the one with the largest last column is ordered above the other 
one. Suppose now that we know the characters of all representations ordered below a representation 
R with Young tableau [ki, . . . , fej]. Consider then the tensor product [ki, . . . , ki_\] (g> [fej]. Both are 
ordered below R and hence their characters are known according to our assumption. The tensor 
product yields [ki, ■ ■ ■ , ki] plus representations ordered below [k±, . . . , ki], and hence we can now 
determine the character of R. Proceeding like this we can systematically compute all characters. 

Not only the characters, but also the trace identities for A r were obtained in p6| 



E ni(-^f-) = ° (k>r + l), (58) 



. m.j\ \ rij 



where the summation is as in equation (pl|). This result was obtained from equation (51) using the 
fact that for A r anti-symmetric tensors of rank larger than the rank of the algebra are trivial. 

To illustrate this let us return to the Weyl formula, equation (|55|). For A\ this yields a very 
simple result for a representation of spin j: 



sinh((2j + l)h) 
sinh(/i) 



Ch^h) = ; V (59) 



Expanding this for the spin-^ representation {j = \) we get 

Ch±(/0 = 2 + h 2 + ±-h* + -U 6 + . . . (60) 
2 12 360 

The spin-^ representation serves as the reference representation in our method. Hence its character 
is 

Tr e F = 2 + \{TrF 2 ) + -^-(TrF 4 ) + ^(TrF 6 ) + . . . (61) 
Using the SU{2) trace identities (||) Tr F 4 = \{TrF 2 ) 2 and 

Tr F 6 = 6(-^-({TrF 2 ) 3 ) + -(TrF 2 ){TrF i ) 
48 8 

1 (TrF 2 ) 3 (62) 
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we arrive at the answer 



Tr e F = 2 + \{2h 2 ) + i(2fc 4 ) + ^(2h*) + ... (63) 

were we substituted F = ha 3 , so that Tr F 2 = 2h 2 (Obviously we could have substituted the 
diagonal form of F directly in equation (^) , but the use of trace identities is far more convenient 
for larger algebras). 

Clearly the equations (^) and ( |63| ) agree, as expected. However, the way the agreement comes 
out is not entirely trivial (although it can easily be derived). Note in particular that the Weyl 
formula is a priori expressed in terms of only r variables, so that all trace identities are already 
built in. On the other hand, in writing down the formal expression Tr e F there is no need to 



18 



specify the number of variables, and indeed the formula is the same for any algebra. The non- 
trivial group structure is thus encapsulated in the trace identities. It is instructive to compare the 
two formulations also for other representations. 

In this case the Weyl formula is superior in elegance and simplicity, although it is somewhat 
more difficult to expand to higher orders due to its denominator. For higher rank groups the Weyl 
formula becomes extremely cumbersome, as explained earlier, while our method does not grow in 
complexity. 

5.4 B r characters 

The basic representations are the anti-symmetric tensors of rank 1 . . . , r — 1 plus the spinor repre- 
sentation. The characters of the anti-symmetric tensors are related to the vector character as in the 
case of A r . The spinor character can be expressed in terms of traces of the vector representation 
by explicit computation. The result is |26]| 



where i?2« are the Bernoulli numbers. 

For algebras of type B r the same trace identity as for A r holds, but with order k > 2r + 1. This 
is true because of the embedding B r C A2 r - All traces of odd order vanish trivially. 

The demonstration that the other characters can be obtained recursively from the tensor product 
rule is similar as for A r , with some complications due to the spinors. We will omit the details (and 
the same holds for C r and D r ). 

5.5 C r characters 

The fundamental representations are the anti-symmetric tensors of rank I = 1 . . . , r with a sym- 
plectic trace removed. The character of the fundamental representation / is equal to the / th anti- 
symmetric tensor power of the vector character minus the (I — 2) th anti-symmetric power of the 
vector character (if I > 2). 

The C r trace identities can be derived using the embedding C r C A^r-x^ which leads to trace 
identities for traces of order k > 2r. Just as for B r , the odd traces vanish. 

5.6 D r characters 

The fundamental representations are the anti-symmetric tensors of rank I = 1 . . . , r — 2 plus the 
two conjugate spinor representations. The anti-symmetric tensor characters are computed as for 
B r , but the spinor characters cannot be expressed completely in terms of traces over the vector 
representation. This is because there exists a symmetric tensor of rank r which never appears 
in traces over the vector representation, namely the Levi-Civita tensor. This tensor is an anti- 
symmetric tensor of rank 2r with vector indices. Combining the 2r vector indices in pairs, with 
each pair labelling an element of the adjoint representation, we can view the Levi-Civita tensor 
also as a symmetric tensor of rank r with adjoint indices. 

Using this new invariant, we can write down the spinor character: 




(64) 




(65) 
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Ch, 



(0,...,0,1) 



(F) 



exp 



E 

_n=l 



(2 



In 



l)B 



--jXr(F) exp 



4n(2n)! 

B2n 



2n Tr F 2n 



oo 



-Tr F 



2)i 



(66) 



where ^yXr(F) is a polynomial of order r in F denned by the leading term in the difference of these 
expressions. It is proportional to the Levi-Civita tensor with indices pairwise contracted with F a . 
The precise definition of the tensor is given in appendix E. 

The trace identities for D r are as those for B r for k > 2r. However, due to the extra fundamental 
trace of order r, there must be an additional trace identity to reduce the number of independent 
ones back to r. Indeed, it turns out the the trace of order 2r can be eliminated using the identity 



e n 



1 



mil 



Xr F ni 



11; 



4(-in± Xr (F)] 



(67) 



This identity was also obtained in [Eq] (the coefficient on the right hand side is incorrect in [26]). 



5.7 G*2 characters 

Exceptional group characters can be computed by expressing them in terms of characters of a 
regular subalgebra. Since the subalgebra has the same rank, one gets polynomials in the same 
number of variables and hence no information is lost. For G2 the only option is the subalgebra A<i- 
We have 

Ch G2)7 = Ch A2 , 3 + Ch A2j 3 + Ch A2il , (68) 

denoting representations by their dimension and omitting the argument F . Since all G2 represen- 
tations are real, the third order invariant of A2 is always cancelled out, and all other odd invariants 
vanish as well. The fourth order invariant can be expressed in terms of second order ones using 
the A2 trace identity. The sixth order invariant of A2 can be expressed in terms of lower ones, but 
this expression involves the third order invariant which doesn't exist in G2- Hence in G2 the sixth 
order invariant is new. After a little algebra we can write the reference character of G2 as 

Ch G2 , 7 = 7 + ^Tr F 2 + ii(Tr F 2 ) 2 + ^Tr F 6 + 1[|(T* F 2 )(Tr F 6 ) - A (Tr F 2 ) 4 ] + . . . (69) 

Here all explicit traces are over the reference representation (0, 1) of dimension 7. The character 
of the other fundamental representation, (1,0) of dimension 14 is easily computed from the anti- 
symmetric tensor product (7 7) a = (7) + (14). Explicitly: 

Ch G2 , 14 = 14 + ^4Tr F 2 + F 2 ) 2 + ^[-26Tr F 6 + ^(Tr F 2 ) 3 ] 

+ ! h ir (T¥ F2)(Tr p6) ~ lr (Tr f2)4] + • • • (70) 

Note that all traces here are over the reference representation. From this expression we read off 
the second and sixth indices of the representation (14): they are 4 and —26 respectively. 

Since the (7) of G2 can be embedded in the vector representation of SO(7), G2 inherits all B3 
trace identities for traces of order 8 and higher. There is an additional trace identity for the fourth 
order trace, which can be read off directly from ChG 2 ,T- 

Tr F 4 = i(Tr F 2 ) 2 (71) 
This exhausts the set of trace identities for G2. 
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5.8 F 4 characters 



The computation is similar to the previous case, now using the sub-algebra B4. We have 

Ch F4j26 = Ch B4i i 6 + ChB 4>9 + Ch B4i i (72) 

The vanishing of the fourth order invariant is not obvious in this case, but follows easily. The sixth 
and eight order polynomials are directly related to those of B4. The tenth order one vanishes again 
by inspection (i.e. the tenth order trace can be expressed in terms of -B4 traces of order 2,6 and 8, 
but not 4), and the twelfth order trace involves the third power of the fourth order polynomial of 
.B4, which did not occur before. It is absorbed in the definition of Tr F 12 , the 12 th order term in 
Chp 4j 26 (up to a factor ^\). To obtain the character of (1,0,0,0) (dimension (52)) we use 

Chi? 4i52 = Ch_B 4 ,36 + C1ib 4i i 6 (73) 

and we express all B4 traces into F4 traces using the definitions introduced when computing [7^. The 
other characters of basic representations can be obtained from the anti-symmetric tensor products 

ChF 4i273 = Chp 4i [ 2 ]*26 — CllF 4i52 (74) 
Chp 4i i274 = Ch F4) [3] !t 26 — ChF 4> 5 2 ChF 4i 26 + Chp 4) 26 (75) 

The last identity involves a little algebra. In the third order anti-symmetric tensor power of 
(26) occurs, in addition to (1274) also the representations (273) and (1053) (with Dynkin labels 
(1,0,0,1)). The former character is known, the latter can be computed using the tensor product 

(26) ® (52). 

Of course the characters of (273) and (1274) can also be computed using the -B4 embedding. 
We have used this check. 

Just as for G 2 one may read off the Dynkin indices from the characters. They are shown in the 
appendix. Furthermore there are trace relations for traces of fourth and tenth order which are read 
off from Ch(26) • By expanding the characters to sufficiently high order one obtains trace identities 
for traces of order 14 and higher. The embedding C -D13 gives trace identities for all traces of 
order 26 and higher, namely precisely those of -D13. We will only present the identities for orders 
lower than that of the maximal Casimir operator. 

5.9 E G characters 

Here we used the sub-algebra Al © A5, and the decompositions 

Ch£ 6i27 = Ch Ali2 Ch A5iB + Ch A5>15 (76) 

Ch B6i78 = Ch4 lj2 + Ch A5i3 5 + Ch Au2 Ch A ^ 20 (77) 

The computation is very similar to the previous cases. We get another basic representation, the 

(27) , by conjugation: 

Ch £ 6 ,27( F ) = Gh E6 , 27 (-F) ■ (78) 

Furthermore the anti-symmetric tensor power of order 2 gives us the representations (0,1,0,0,0,0) 
(351) and (0,0,0,1,0,0) (351) and the order three anti-symmetric power yields precisely the repre- 
sentation (0,0,1,0,0,0) (2925). 

The indices and trace identities for orders up to 12 are listed in the appendix. 
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5.10 E 7 characters 

Here we used the sub-algebra Aj, and the decompositions or anti-symmetric tensor products 



ChE 7i56 = Ch A7>28 + Ch Ar ^ (79) 

Ch S7il33 = Ch A7i70 + Ch ArM (80) 

Ch^ 7 ,9i2 = Clu 7 ,420 + Ch A7 + Clu 7 ,36 + Ch^gg (81) 

Ch^ 7i i539 = Ch S7i [2]*56 — Ch^ 7) i (82) 

Chg 7 ,8645 = Ch S7) [ 2 ]*133 ~ Ch£; 7> i3 3 (83) 

Che 7i 27664 = Ch_E 7i [ 3 ]^56 — Ch£ 7)56 (84) 

ChE 7 ,365750 = Ch^ 7j [4]*5 6 — CriE 7i i539 — CllE 7) i (85) 

5.11 E 8 characters 

Here we used the sub-algebra D$, and 

Ch£ g ,248 = Ch£) gj i28 + CllD g ,120 (86) 

ChE gj 3875 = Ch/j gj i920 + Cri£) 8i i820 + Cho 8i i35 (87) 
Ch£; g , 147250 = Ch/j gi 60060 + Ch/j gi 56320 + Ch/j gj i5360 

+Ch£) gi 7020 + Cho g) 6435 + Ch/3 gj ig20 + CliD gj i35 (88) 

Ch£; 8i 30380 = Ch Sgi [ 2 ]*248 ~~ Chs g ,248 (89) 

Ch^g, 2450240 = Ch Egi [31^248 — (Ch Sgi2 48) 2 + Ch Egi2 48 (90) 

ChE gi66 96000 = Ch Sg) [ 2 ]*3875 — Ch^ gj38 7 5 (ChE g ,248 — 1) + Chs g) i47250 (91) 

ChE gi 146325270 = Ch Eg) [4]* 2 48 — (Ch Eg! [ 2 ]*248 — CllE g) 248) (ChE g)2 48 ~ 1) (92) 



ChE g)689 9079264 — Ch^ 8) [ 5 ]*248 ~~ Chs g ,248(Ch £ ; gj [3^248 — 2Ch Bgj [ 2 ]*248 + Ch^ g) 248 ~ 1) (93) 

In the first three lines all D$ spinor representations must be from the same conjugacy class, which 
is fixed by the decomposition one chooses for the (248). Since the choice one makes for the class 
is irrelevant, there is no need for a label to distinguish conjugate spinors. The representation 
(6435) is an (anti-)selfdual tensor. It belongs to the trivial conjugacy class, but it does carry a 
non-trivial chirality. If for the representation denoted (128) we choose the one with Dynkin labels 
(0, 0, 0, 0, 0, 0, 0, 1), then the correct set of Dynkin labels for the 6345 is (0, 0, 0, 0, 0, 0, 2, 0). 

5.12 Normalization of indices 

The normalization of the symmetric tensors is fixed by fixing a normalization for the indices. We 
will do this in such a way that they are always integers, as the word "index" suggests. For the second 
index there is a natural normalization in terms of the Atiyah-Singer index theorem for instantons 
on 53. For any representation of any algebra we can choose the second index equal to the net 
number of zero modes of a Weyl fermion in that representation in an instanton field of minimal 
non-trivial topological charge (where "net" means the difference between the two chiralities) . Then 
the second index is equal to 1 for the reference representations of A r and C r , 2 for those of B r and 
C r , and 2,6,6,12,60 respectively for the reference representations of G2, Eq, £7 and Eg. For the 
adjoint representation the second index is always equal to twice the dual Coxeter number g listed 
in table [|. This choice corresponds to setting 77 = 2 in (|27]). This value of 77 was used in the last 
column of table ||. 
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For the higher indices there is a similar topological interpretation in terms of gauge bundles 
on higher dimensional manifolds, but we will not explore that here in detail. One may however 
follow the spirit of such an interpretation and define all higher indices in such a way that they are 
integers. This is automatically true if they are integers for the basic representations, because the 
characters of all representations are polynomials with integer coefficients in terms of the characters 
of the basic representations. Furthermore, within the set of basic representations the ones obtained 
by means of anti-symmetric tensor products of a given representation R, have indices that are an 
integer multiple of those of R. Then only the spinor representations of SO(N), one representation 
of F4, Eq and E-j and two of Eg require special attention. 

For the reference representations we choose all higher indices equal to 1, except when a larger 
integer is required to make all indices integral. The 2n th index of a spinor representation of SO(N) 
follows directly from the character equations (^) , ( |65| ) and (|66|) : 

d, m( 5) ' 22 "- n 1 '^ , (94) 

where dim(S') is the dimension of a spinor representation, i.e. dim(S') = 2 r for algebras of type 
B r (SO(2r + 1) and dim(S) = 2 r " 1 for type D r (SO(2r)). This assumes that the (2n) th index of the 
vector representation is set to 1. By inspection, this expression is an integer except for the fourth 
index for SO(N),N < 8 and the eighth index for N < 10. We have checked that the spinor index 
is an integer in all other cases for 2n < 100. Table ^| gives the indices for the SO(N) vector and 
spinor representations, according to our normalization (the index I r is not listed here; its value 
is for the vector representation and chosen ±1 for the two fundamental spinors of SO(N), N 
even). Note that for N = 3, ... ,6 the spinor representation, and not the vector is the reference 
representation, which automatically leads to the entries in the table. 
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h 


h 


h 


h 


ho 


h2 


hi 


3 


4,1 














4 


4,1 














5 


2,1 


-4,1 












6 


2,1 


-4,1 












7 


2,2 


2,-1 


1,1 










8 


2,2 


2,-1 


1,1 










9 


2,4 


1,-1 


1,2 


2,-17 








10 


2,4 


1,-1 


1,2 


2,-17 








11 


2,8 


1,-2 


1,4 


1,-17 


1,124 






12 


2,8 


1,-2 


1,4 


1,-17 


1,124 






13 


2,16 


1,-2 


1,8 


1,-34 


1,248 


1,-2764 




14 


2,16 


1,-4 


1,8 


1,-34 


1,248 


1,-2764 




15 


2,32 


1,-4 


1,16 


1,-68 


1,496 


1,-5528 


1,87376 


16 


2,32 


1,-8 


1,16 


1,-68 


1,496 


1,-5528 


1,87376 



Table 4: Indices for the SO(N) vector and spinor representations. The chiral index I r is not 
listed here; we choose it equal to +1 for the fundamental spinor (0, . . . , 1, 0) and equal to —1 for 
(0,...,0,1). 

For the exceptional algebras the two highest indices of E7 and the three highest ones of Eg 
come out fractional unless we choose a different normalization for the reference representation. 

Our preferred index normalization for the reference representations is summarized in the last 
column of table ||. For SO(2N) the last entry indicates the index 1^ of the spinor representation 
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(0, . . . , 1,0). Only the second index is affected by the choice of ij. For all higher order traces we 
fix the normalization of the index, and then the ^-dependence goes into the normalization of the 
symmetric tensor. 

This also fixes the normalization of all symmetric tensors as 

Str T r ai . . . T?" = I n (r)d ai - a " . (95) 

This defines the tensors given a choice of generators in the reference representation. 
In Appendix D we present some results for contractions of these tensors. 



6 Reduction identities 

At this stage we have terms that contain combinations of the invariants d r and the structure 
constant /. Our task is now to eliminate / from the terms as much as possible, and reduce the 
total number of invariants in the final answer. 

Unlike the results obtained so far, for these reductions we cannot give a general algorithm. In 
fact, we do not even know what the desirable outcome is, since we are not aware of a mathematical 
theorem that gives us a basic set of invariants in terms of which all others can be expressed. At 
any given order we can derive large numbers of identities among the various invariants, but there 
will be new relations at every order. In practice this is not a major problem. First of all we have 
obtained results relevant for vacuum bubble Feynman diagrams of up to nine loops, and secondly 
the number of invariants we are left with is small, although possibly not minimal. 

The most useful identities for doing this are the Jacobi identities. Thus far we have seen two 
of them in equation (fl3l). The second equation there can be generalized into 

cyclic permutations of i\---i n 

for all representations R, including the adjoint representation. In fact, there is an advantage 
to postponing the replacement of the adjoint symmetric tensors by reference tensors, since they 
satisfy additional identities. Furthermore some identities produce new tensors d A - For this reason 
we present the results in terms of dn and d A rather than d r . 

The first identity that can be derived from this is one for invariants with three indices: 

d< M faj c f blc = \ CA df (97) 

It is actually the simplest identity in a class of triangle reductions that involve one or two invariant 
tensors d^. We have also 

jixjx—jnkx jiijx— j n kl ekxk%%z 
Rx R2 J 

j>xix-"jnkx ffiS— imjx-'jnk2 fkxkzh 

Ri R2 J 

These identities are very powerful when a large number of invariants is involved. 

For showing further reduction identities we will use a special notation which corresponds closely 
to a notation that can be used inside a computer program. We will represent a (symmetric) invariant 
by a product of vectors: 

d%- in = p^---p{ n (100) 



1 1 



Na n + l Rl 



jjjx—On+2 ^jx—jn+2 Hii2h 
Ri R2 J 



n + l 



-d 



OX- -jn + lk ^13 ■ --imjl ■ ■ -jn + 1 jkill2 



Rx 



d 3 ' 
a R 2 



(98) 
(99) 
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The lower index on the vector refers to the particular invariant. In this notation we have no 
problems with the symmetric property of the invariants. Of course we are not implying that each 
invariant can be mathematically written this way. It is just notation. 

Additionally we will use Schoonschip notation on contracted indices. That is: if the index of a 
vector is contracted with an index of a tensor we put the vector in the place of this index. Hence 



f PlP2i f PlP2i (pi- P 2) 



(101) 



Furthermore we can add a weight to each formula. This weight is basically the number of vertices 
(assuming that all vertices are three-point vertices) in the diagram before we started the elimination 
procedures. For our current algorithms the weight is the total number of vectors pi plus the number 
of structure constants /. Hence the weight of the above formula is 2n + 6. We will present all 
reduction identities that are relevant for weights up to 12. This corresponds to 7-loop vacuum 
bubbles or 6-loop propagator diagrams. The derivation of all these identities involves the use of 
the generalized identity of equation @). 

jP\Pli\ JPlP2«2 JPliliz /P 1 * 2 * 3 (p 1 -p 2 ) n 
JPlP2«i jPxhil jVli-ziz fPlPzh ^ .p 2 \ n 



1 



n + 1 
1 

n + 1 



Ca (P1-P2 



>n+2 



dP A 



1P1P1P2 



(Pi -P2 



.n+1 



6(ra+l)(n + 2) 



Ca (prP2 



,n+3 



1 



n + 



-dPf P2P2 {p v p 2 ) 



n+1 



ypiP2U jpiP2«2 ypiP2«3 f iliii3 (j)i-p2) r 

jPlP2h JPlP2h JPlP2«2 f PlV2i2 (p l -p 2 ) r 



f PlP2i f PlP2i (pi-P3)(p2-PsT 
f PlP2i f PlP3l (pi-P2)( P 2-P3T 
f PlP2i f PlP3l (pi-P2)(p2-P3)(pi-P3T 

f plP2m f plP2P3 (pi-P2)(prP3) 







9 

_ n - JP\P\PZP2( m _ \n+2 

(^+2)1^ (P1 ' P2) 



3 (n + 3)! 



Ca {pi-P2) 2 {pi-P3){P2-P3) n 
Ca (pvP2) 2 (pi-P3)(P2-P3) n 



C A 



2(n + l) 
1 



{prP2) 2 {P2-P3){prP3 



.n+1 



^ Ca (prP2) 2 (prP3) 2 (P2-P3) 



(102) 
(103) 

(104) 
(105) 

(106) 

(107) 
(108) 
(109) 

(110) 



The equations with an odd number of f's that are relevant are all zero. This has been explicitly 
shown, but for many of them one can see this already on the basis of symmetry principles. 

We are not going to present the identities that would be needed for diagrams of weight 14 or 
16. There would be too many of them and moreover, this is not how we have constructed the 
computer program. In the program we have found a way to apply equation (|96| ) recursively in such 
a way that it reduces all combinations with the exception of one (at weight 14). The program has 
to guess at what combination of invariants and /'s to take for the application of the formula and 
we let it guess several times. In the end this covers all cases except for the one that we cannot do 
by these methods anyway. 

The first object that causes some real problems because the above algorithms are not sufficient 
to handle them, occurs at weight 14. This object can either be written as 

f PlP2P3 f PiP2P3 (pvP2)(pvP3)(p2-P3) (HI) 
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or with some rewriting (and omitting trivial terms that are a byproduct of the rewriting): 



f PiP2i f P1P2i (pi-p3) 2 (p2-P3? (112) 

If at least two of the three invariants are in the adjoint representation this object can be reduced 
with the same technique that we use below to simplify some combinations of invariants only (see 
equation ( |119|) ). 



6.1 Combinations of invariants 

Here we will consider combinations of invariants only. The easiest combinations are full contractions 
between two invariants as in 

d X' ind %' in ( 113 ) 

Unfortunately, when the weight of the diagrams increases, the complexity of the combinations 
increases correspondingly. In some cases one can make reductions. For instance: 

4; 1 '• = 6 jk ,r^- '■ d^y ^ (114) 

and hence 

m h—in.J) h—in^a ji—j m J> ji—j m 1 J.oii-'-in jioii—in jjoji—jm jjoji—jm (115) 

Ri R2 R3 R4 Na ^ 1 ^ 2 ^ 3 ^ 4 

but for objects of the type 

d%i 2i3i4 d%i 2h d%l 4i5 (116) 

or 

^ili2i3i4«5^*l*2«3«6 ^425*6 (117) 

there does not seem to be a general simplification of this type. In the case of all invariants belonging 
to the adjoint representation we can still do things as we see in the next formula: 

d ^cdef d abcdef_^ d abcd d cdef d ejab + ^ = Q (ug) 

The derivation of this formula is a matter of evaluating a circular ladder with 6 rungs in two 
different ways. In the first way one sees it as two loops with 6 vertices, and in the second way as 
three loops with 4 vertices. The algorithms of the previous sections are then sufficient to obtain 
this formula. Note however that such derivations usually need the use of a computer program: the 
intermediate stages can contain large numbers of terms. A similar technique can be used for the 
object in equation (|112|) . We look at a circular ladder with 7 rungs. If we see this as three loops 
with 4 vertices (and two /'s left) we get the form of the equation, and if we see it as two loops 
with 7 vertices, we get a representation involving two dA invariants with 6 indices (the ones with 
7 indices are zero for the adjoint representation). The result is (after a rather lengthy calculation, 
applying equation ( |11§| ) and normalizing): 

o 

jabcd icdef ,efgh rani fbhi , sy3 rabcd jabcd 
R A A J J + ~^ U A a R a A 

-^C A dt d dT f dT b + l d T defd R bCg d d A fg = (H9) 

with R replaced by A. Because a similar exchange of the order of evaluation in one of the examples 
below gives the same equation with a slightly extended generality we have already presented this 
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more general form with the representation R. One can also derive this equation from equation 
(|14§| ), but the derivation of that formula uses basically the same technique. 

Considering that the algorithms we have presented can reduce all combinations of invariants 
and /'s, with the exception of the above combination with at least two invariants not in the adjoint 
representation, up to weight 14 we have only a limited number of topologies (contractions between 
invariants) left. We will show them graphically, omitting the very trivial ones that can be reduced 



with equation |114| and objects of the type d R 1 " n . The elements of this pictorial language are 



Y 



J«l«2«3 

a R 



J«1*2«3M 

a R 



J«l«2*3«4 

a A 



J«3*2«l 



(120) 



in which we assume the indices of / to run counterclockwise in the diagram. For weight 6 we have 



For weight 8 there is also only one topology: 

^44(Pl,P2) : 

For weight 10 there are two topologies: 

^433(Pl,P2,P3) 

For weight 12 we have 5 topologies: 

dm(pi,P2) 

de33(Pl,P2,P3) 
d543(Pl,P2,P3) 

du4(pi,P2,P3) 

d3333(Pl,P2,P3,Pi) 



<X) 3 jg) 
1 2 



4 (X ) 

1 2 



(X) 5 (X ) 
1 2 




1 2 

2 13 




(121) 

(122) 

(123) 
(124) 

(125) 
(126) 

(127) 
(128) 
(129) 
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Finally for weight 14 we have 9 topologies: 

d743{Pl,P2,P3) 
d653(Pl,P2,P3) 

d644(Pl,P2,P3) 



7 

1 2 

2 13 



(130) 
(131) 

(132) 
(133) 



d554(pi,P2,P3) = 

d5333{Pl,P2,P3,P4) = 

du33a(pi,P2,P3,P4) = 

d4433b{Pl,P2,P3,P4) = 

du33c{Pl,P2,P3,P4) = 






(134) 
(135) 

(136) 
(137) 
(138) 



For weight 16 there are more than 20 non-trivial topologies. We will not show them here. 

In terms of these diagrams of invariants we can compose a few extra equations that can be very 
useful: 



(Cr - -C A ) I 2 {R) 5^ 
6 

(c R -\c A )^ 

(C R -\c A )(^z +^h{R) 
(C R --C A )^ +- 



<8b^ = -\c A -ngng^ - ^h{Ri)C 2 A 

2 l""-- 8 2 I"""- 40 2^ 

= -\ Ca + ivtr?^ 



i 



N Rl 2 1 1 
43 1 



120 N A 21 



23 
1440 



2 1 



49 
60 



-©- 




(139) 
(140) 
(141) 
(142) 
(143) 
(144) 

(145) 



We assume symmetrization over the external legs in the terms of the right hand sides. 
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One of the spin-offs of equation ( |144| ) is 
2 13 o 




^/ 2 (i?i)Ci(8B8) + -i- 

40 2 3 N Rl 2 



1 1 3 

7 fc 



(146) 



This allows us to eliminate the topology 0^33 completely. For groups for which d lJ exists and 14(A) 
is not zero we can use the techniques of the next section to also eliminate the topologies 0(433 , 0(44330 
and g?4433 C . Similarly equation (|145| ) gives 



(X) 4 C*5=i%) 
3 12 




43 1 in ^ n 

— — <8B8) mm 

120 iVA 2 1 3 



(147) 



and therefore also the topology 0(743 can be eliminated. 
An equation that is also rather interesting is 

:C 3 




(148) 

It can be used to derive equation (|119| ), but it can also be useful for bigger diagrams. It is derived 
by writing a diagram with 10 /'s into loops in two different ways. After that the application of the 
reduction algorithms and some rewriting leads to this formula. 



7 Representation independent invariants 

The invariants that we have used thus far were only symmetrized traces. It is possible to define a 
new set of invariants that is not only symmetric, but also orthogonal (see equation (21)). Because 
the invariant with two indices is proportional to the Kronecker delta in the adjoint space, this 
means for instance that these invariants have a zero trace (a contraction of any two indices gives 
zero). But also contractions with all the indices of invariants with fewer indices than the invariant 
under study should give zero. For most algebras we can define (the exceptions are certain SO (AN) 
algebras with two independent tensors of order 2 N) : 



d% i2isk = h(R) d^ 2 ^ + h,2(R) (5 lll2 5 lsli + 5 llls 5 l2iA + 5 ni4 5 i2i:i )/3 (150) 
d ni2i 3 i4i 5 = j 5 ( E ) rfiifaisuis + I^ 2 (R) (rfiiwajuw + . . .)/ 10 ( 151 ) 

d hi2huisi 6 = d iii 2 w 5 *6 + I^R) (^sU^ie + . . .) / 15 

+I 3>3 (R) (df^df 5 ^ + ■ • -)/10 + I 2 ,2,2(R) (5 Hi2 5 i3h 5 i5ie + • • -)/15 (152) 
^Maiskfewr = j^isuiswr + I 5 2 (R) ^ivtiaUh giair + . . .) / 2 1 

+h,s(R) (df^df^ 7 + • • -)/35 + h,2,i(R) (d^d 1 ^^ 7 + • • -)/105(153) 

and higher ones are defined analogously. The composite constants like l2,2(R) can now be derived 
by the orthogonality conditions. When we multiply equation ( |150| ) by 5 laU we obtain 

Na + 2 ] I 2 ,2(R)5 hi2 = df 2i:ii3 
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Similarly we can derive 



= (C R -C A /6)4 h 
= (C R -C A /6)h(R)5 ili2 
^^h,2(R) = (h(R)^ R -\h{A))h{R) (154) 



h,2(R) = j^^{C R -\c A )h{R) (155) 



For the next level of constants we get coupled equations. The easiest way to derive them is with 
the use of the equations (|142|) and (|146|). If we contract the first equation once with efi* 2 * 3 * 4 and 



once with S ni2 5 l3U and substitute equation (152) in all three of the equations that we obtain this 




way we have enough to come to a solution. Unfortunately this solution is not very elegant. In 
many cases one can make simplifications. For instance one can write (using equation ( |139| )) 

= \c\ ®=® 

6 

= h(A) + \h,2{A) ®=® (156) 

which then reduces to 

We notice that whenever Ia(A) = or cfi* 2 * 3 * 4 = we must have that either = or 

N A = 8 (assuming that Ca is never zero). The last is indeed the case for SU(3). Hence a general 
application of this relation is rather dangerous. 

The equations for 1^2 (R), ^3,3 (R) and I2,2,2(R) can be derived by taking the definition of 
the decomposition in equation Q152Q and multiplying either by 5 Hl2 5 l3U 5 t516 or rf^ 21 * 1 * or 
<^l i2h dtf 5i \ With the use of the equations O and @) we obtain: 



= 0=0 &hfl{R)) 
5 

+N A (N A + 2){^-^I 2X2 {R) - \{C R - \c A )h.2{R) - ±-h,2{A)h{R)) (158) 
15 3 3 90 

(\hj{R)) 




+ ©Sr) (^±^J 4i2 (fl) - (C R - \c A )h(R) - ±-h(A)h(R)) (159) 
15 3 30 

[\h,2{R) + \c A h(R)) 




^^n^hs(R)-^iHR)?) 

+ 0=^r) &2,2,2(R) + \c A I 2 ,2{R) + - A C\h(R)) (160) 
5 4 4 

Because for each group either some of the objects in these equations are zero, or there are simplifi- 
cations, it does not seem wise to solve this system in this form. For all groups with the exception 
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of SU(N) we have that d 1 ^ 213 = 0, and hence the system reduces to two equations in /42CR) and 
^2 2,2 (R)- For SU(3) we have that d % ^ 2l3%i = and again we have a simpler system. For the other 
SU(N) groups we can apply equation ( |157| ) improving the solutions somewhat. Because of the 
singularity (zero divided by zero) of the solutions for all groups but SU(N),N > 3, we do not 
present the general solution here. It serves no purpose because we will not use them. Along the 
same lines one can derive the equations for Iq^R), ^4,3 (R) and h t 2,2(R)- 

The expressions at rank 6 and 7 become rather complicated due to the fact that the tensors in the 
right hand side of the equations ( [L52D and (|l53j ) are not orthogonal. Hence the various orthogonality 
relations mix and this gives the complicated result. It should be clear that invariants with a higher 
rank will give even more complicated relations. The exception is the adjoint representation. For 
this representation all invariants with an odd number of indices are zero. Hence /3(A) = I5 (A) = 

/ 3 , 2 (A) = ... = 0. 

We can use the above relations for the reduction of some contractions of invariants. We have 
seen in the previous section that for the adjoint representation more things are possible than for 
the other representations. However we can rewrite the invariants only to invariants of the adjoint 
representation when the corresponding cPs can indeed be expressed as such. This means that 
the corresponding /(A) should never be zero. Unfortunately this excludes many contractions of 
invariants. Hence we do not see many benefits here. 



7.1 Orthogonal versus Reference tensors 



Here we will compare the computation of a symmetrized trace in two ways, using the orthogonal 
basis (satisfying (|2l|)) and using characters, with tensors defined for a reference representation. The 
latter will be referred to as "reference tensors" . We consider fourth order traces in SU (N) . In this 
case an explicit expression exists for any representation 



Str TrTrTrTr 
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(161) 



(162) 



Here "d" without subscript denotes the orthogonal tensor. The reference representation is the 
vector, and the reference fourth rank tensor is by definition equal to ( |l6l| ) with "i?" equal to the 
reference representation: 
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jabed 
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(163) 



(Note that ^2,2 is the same in both bases). This allows us to express (161) in terms of reference 
rather than orthogonal tensors: 

/ 4 (RK bcd 



Str T%T b R T c R T* 



N A 



h{Rf 



Na 
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"2,2 



(164) 



Consider now for example the anti-symmetric tensor representation. Its character, expanded to 
fourth order is 



Ch [2] (F) 



1 



N(N - 1) + -(N - 2)Tr F 2 + -(N - 4)Tr F 3 
2 6 



+— (N - 8)Tr F 4 + -Tr F 2 Tr F 2 + . . . 
24 v ' 8 



(165) 
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where all traces are over the vector representation. From the fourth order terms we deduce, by 
differentiating with respect to F a , . . . , F b : 

Tr r^r^T^T^] = (N - 8)df cd + 3d^ d (166) 

We may now verify this using ( |l64f) . Although the coefficient of the second term does not look 
very encouraging, substituting hiR) = N — 2, Na = N 2 — 1 and Nr = ^N(N — 1), it does indeed 
produce the coefficient 3 in (|166|). 



This illustrates several points. Trace formulas in terms of orthogonal tensors such as (161) have 
a simpler form than those in terms of reference tensors, if one tries to write down expressions for 
arbitrary representations R. However, expressions such as ( |166| ) can be written down fairly easily 
for any representation although not (easily) in closed form. Furthermore they can be extended to 
arbitrary order in a straightforward way while this rapidly becomes extremely difficult for ( |161| ) or 

dH). 

Note that the indices (the coefficients of fundamental tensors) are basis independent (apart from 
normalizations), whereas the sub-indices (coefficients of combinations of fundamental tensors) are 
not. In the orthogonal basis it is not hard to see that all sub- indices can in fact be expressed in 
terms of indices, so that they do not constitute an additional set of variables. In any reference 
basis the same is then true, since it can be related to an orthogonal basis, but the expression are 



(even) more complicated, as in (164). For orders larger than six, expression of sub-indices in terms 
of indices are not available and hard to obtain. 

In our method the sub-indices are essentially treated as additional variables, which can be 
computed for any representation as easily as the indices themselves. This allows the computation 
of a trace for any representation, which was our goal. The result is a combination of symmetric 



fundamental tensors with explicit numerical coefficients, as in (166), or an expression involving 



both indices and sub-indices. Unfortunately it is much harder to present the result in minimal 
form, with all representation dependence encapsulated in the indices. 

A Indices and Trace identities for Exceptional Algebras 

In this appendix we summarize our results on traces for exceptional algebras. All the indices for the 
lowest dimensional representations are given, including all basic representations. Trace identities 
are given for all traces of order less than the dual Coxeter number g. 

For the classical algebras A. . . D the trace identities were already given in chapter 6, and index 
formulas for some representations are given in ||. 

The indices provide only a small part of the information contained in the full characters, but it is 
impractical to present the latter in printed form. We do have an efficient procedure to generate the 
characters of any representation of any simple Lie-algebra to any desired order. This procedure uses 
a combination of Kac [^lj (to compute tensor products) and FORM [^] (to multiply, add and subtract 
characters according to these tensor products), and is available via http: //norma. nikhcf.nl/^t5<: . 



Obviously this then also provides all indices for algebras and representations not listed in this 
appendix. 

A.l Indices and trace identities for G 2 

Trace identity in the representation (7): 

Tr F 4 = i(Tr F 2 ) 2 (167) 
The indices of the lowest-dimensional representations are shown in table 5. 
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Rep. 


Dimension 


2 


h 


(0,1) 


7 


1 


1 


(1,0) 


14 


4 


-26 


(0,2) 


27 


9 


39 


(1,1) 


64 


32 


-208 


(0,3) 


77 


44 


494 


(2,0) 


77 


55 


-1235 


(0,4) 


182 


156 


3666 


(1,2) 


189 


144 


-456 


(3,0) 


273 


351 


-20709 


(2,1) 


286 


286 


-7904 


(0,5) 


378 


450 


19500 


(1,3) 


448 


480 


2640 


(0,6) 


714 


1122 


82212 


(2,2) 


729 


972 


-27378 


(4,0) 


748 


1496 


-193324 


(3,1) 


896 


1472 


-109408 



Table 5: Indices for G<i- 



A. 2 Indices and trace identities for F A 

Trace identities in the representation (26): 

Tr F 4 = 3(^Tr F 2 ) 2 (168) 

Tr F 10 = jj^Tr F 2 )(Tr F 8 ) - ^Tr F 2 ) 2 (Tr F e ) + ^r(^Tr F 2 f (169) 
The indices are listed in table 6. 



Rep. 


Dimension 


h 


h 


/ 8 


hi 


(0,0,0,1) 


26 


1 


1 


1 


1 


(1,0,0,0) 


52 


3 


-7 


17 


-63 


(0,0,1,0) 


273 


21 


1 


-119 


-1959 


(0,0,0,2) 


324 


27 


57 


153 


2073 


(1,0,0,1) 


1053 


108 


-132 


612 


372 


(2,0,0,0) 


1053 


135 


-645 


2907 


-134373 


(0,1,0,0) 


1274 


147 


-133 


-1309 


125811 



Table 6: Indices for F4. 



A. 3 Indices and trace identities for E 6 

Trace identities in the representation (56): 

Tr F 4 = 12(— Tr F 2 ) 2 
v 12 ; 

Tr F 7 = ^(Tr F 5 )(^Tr F 2 ) 



(170) 
(171) 



Tr F 10 = -(Tr F 8 )(— Tr F 2 ) - 7(Tr F 6 )(— Tr F 2 ) 2 + -^(Tr F 5 ) 2 + 42(— Tr F 2 

^ .1. ^ .1_ L. W ) _1_ _ 

Tr F 11 = -(Tr F 6 )(Tr F 5 ) + — (Tr F 9 )(-Tr F 2 ) - -(Tr F 5 )(-Tr F 2 ) 3 
The indices are in table 7. 



Rep. 


Dimension 


h 

6 


h 


h 


h 


^9 


Il2 


(1,0,0,0,0,0) 

\ 7 7 7 7 7 / 


27 


1 


1 


1 


1 


1 


1 


(0,0,0,0,1,0) 

\ 7 7 7 7 7 / 


27 


1 


-1 


1 


1 


-1 


1 


(0,0,0,0,0,1) 

\ 7 7 7 7 7 / 


78 


4 





-6 


18 





-62 


(0,1,0,0,0,0) 

\ " 7 7 7 7 / 


351 


25 


11 


-5 


-101 


-229 


-2021 


(0,0,0,1,0,0) 

\ 7 7 " 7 7/ 


351 


25 


-11 


-5 


-101 


229 


-2021 


(0,0,0,0,2,0) 


351 


28 


-44 


58 


154 


-284 


2074 


(2,0,0,0,0,0) 


351 


28 


44 


58 


154 


284 


2074 


(1,0,0,0,1,0) 


650 


50 





60 


36 





116 


(0,0,0,0,1,1) 


1728 


160 


-88 


-80 


664 


152 


424 


(1,0,0,0,0,1) 


1728 


160 


88 


-80 


664 


-152 


424 


(0,0,0,0,0,2) 


2430 


270 





-720 


3672 





-131928 


(0,0,1,0,0,0) 


2925 


300 





-270 


-918 





122202 



Table 7: Indices for F 6 . 



A. 4 Indices and trace identities for E 7 

Trace identities in the representation (56): 



Tr F 4 = 24(— Tr F 2 ) 2 



16 = _£ju/_,_i_ r F 2 )(TrF 6 )(T r F 8 ) 
5220 v 24 A A ' 
2360 ^ 2 14 

319 v 24 A ' 
61607, 1 m _ 9 . 9/rT , ^ fi ,o 

+ 23490(24 Tri? )( Tri? ) 
_63700 p2 2 12 

2871 v 24 ; v ; 
21164 1 2 3 10 

783 v 24 1 y 1 



7397 p2 4 p8 

522 v 24 ' v ' 
72254, 1 m ~^ /rT , 



783 '24 



Tr F 2 f(Tr F b 



222898 . 1 m o x8 

+ (— Tr F 2 8 

261 v 24 7 

1 3 

+ ^(Tr F 6 )(Tr F 10 ) 



The indices are in table 8. 
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Rep. 


Dim. 


7 1 

12 


h 


h 


Iw 


h2 


hi 


/l8 


(0,0,0,0,0,1,0) 


56 


1 


l 


1 


1 


1 


29 


1229 


(1,0,0,0,0,0,0) 


133 


3 


-2 


10 


-2 


-30 


542 


-111658 


(0,0,0,0,0,0,1) 


912 


30 


-10 


-82 


230 


-2082 


-39170 


96018190 


(0,0,0,0,0,2,0) 


1463 


55 


90 


174 


570 


2134 


238650 


161267970 


(0,0,0,0,1,0,0) 


1539 


54 


24 


-72 


-456 


-1992 


-235944 


-161018664 


(1,0,0,0,0,1,0) 


6480 


270 


30 


774 


-210 


534 


73350 


-102108810 


(2,0,0,0,0,0,0) 


7371 


351 


-354 


2682 


-834 


-63438 


4748094 


-14489069226 


(0,1,0,0,0,0,0) 


8645 


390 


-200 


40 


760 


57480 


-4368520 


14620498520 


(0,0,0,0,0,3,0) 


24320 


1440 


3600 


10176 


50160 


292896 


59512080 


167838228720 


(0,0,0,1,0,0,0) 


27664 


1430 


-10 


-3442 


-7450 


63998 


32976190 


149694252430 


(0,0,0,0,0,1,1) 


40755 


2145 


530 


-3658 


13490 


-171138 


2436850 


-9081228710 


(0,0,0,0,1,1,0) 


51072 


2832 


2872 


256 


-16568 


-172464 


-46178632 


-158703316792 


(1,0,0,0,0,0,1) 


86184 


4995 


-3165 


963 


36195 


-366717 


-37725705 


-137019575865 


(1,0,0,0,0,2,0) 


150822 


9450 


8400 


41328 


59280 


410928 


30093840 


30366263760 


(1,0,0,0,1,0,0) 


152152 


9152 


-328 


9320 


-78088 


-197560 


-28617992 


-27126731432 


(3,0,0,0,0,0,0) 


238602 


17940 


-26380 


271676 


-116620 


-13615284 


1492228660 


-16354668799100 


(0,0,0,0,0,0,2) 


253935 


17820 


-9000 


-94824 


404280 


-6024744 


-323856360 


12685209865560 


(0,0,0,0,0,4,0) 


293930 


24310 


88400 


329888 


2153840 


17066368 


4880546320 


30243257914480 


(2,0,0,0,0,1,0) 


320112 


21762 


-9318 


155826 


-75318 


-3178974 


303759378 


-674257133022 


(0,1,0,0,0,1,0) 


362880 


23760 


600 


12672 


22440 


3531792 


-239671080 


806089955880 


(0,0,1,0,0,0,0) 


365750 


24750 


-9000 


-63240 


79800 


2601720 


278208600 


-12473996293800 



Table 8: Indices for E 7 . 



A. 5 Indices and trace identities for E 8 

Trace identities in the representation (248): 



Tr F 4 = 36(^Tr F 2 ) 2 

Tr F % = 30(^Tr F 2 ) 3 

Tr F w = -(Tr F 8 )(-Tr F 2 ) - — (-Tr F 2 ) 5 
4 v A 60 ' 4 v 60 ' 



(176) 

(177) 
(178) 



Tr F 16 = 



143 /m ^ Sx a 

+ Tr F 8 ) 

1920 v 1 

— ^^(Tr F 12 )(—Tt F 2 ) 2 
20 v n 60 1 

+ ^( Tr F 14 )(-Tr F 2 ) 
5 v A 60 ' 

147147, 1 m 

-^(eo Tr 



(179) 



Tr F 22 = 



+ 



-^^(Tr F 8 )(Tr F 12 )(-Tr F 2 ) 
39360 v n n 60 1 

323 



(Tr F s )(Tr F 14 ) 



1920 



35 



3193461271 /rTn n8u l ~. 7 

+ -f679360- (Ti 60 
361151791 J_ 

50380800 v ' v 60 ; 

- 298544701 (Tr F")(l Tr F>)= 
524800 v A 60 7 

3438981 /m , 4w 1 m „ 2 , 4 

+ Tr F 14 — Tr F ) 

16400 v A 60 7 

-i^Z ( TrF 18 )(l T rF 2 ) 2 
6560 V A 60 7 

164475086139 1 211 
3358720 l 60 J 

1888509563603 8 12 1 2 3 

______ (Tr F )(Tr F )(-Tr F ) 

13299145169 , m o., m „ 14w 1 m „,x2 
- 1347978240 (Tr 60 

400077914058102709 /rTn 1 m ». q 

Tr F 8 )( — Tr F ) 

7074189803520 v A 60 ; 

10415984073858649 , m o, 2 . 1 m 

Tr F 8 ) 2 (— Tr F ) 

35370949017600 v ; v 60 ' 

15835512121 J_ 
221838704640 v ; v 60 ' 

+ ^_( T rF 12 )(Tr F 14 ) 
42336 v A ' 

96960058189033 /rTn , 2w 1 r 
_i /fry rr^v r 

5554483200 1 A 60 

_3481855 J_ 

6318648 v ; v 60 ' 

28373592046607 /rTn , 4w 1 

4386278400 1 A 60 ; 

1446721465417 , m ,«,., 1 , 

+ (Tr F 18 )(— ' 

1973825280 v A 60 

_ 1354577341 _1_ 

7441008 v A 60 

929825 . m o 4w 1 m o, 

+ Tr F 24 )( — Tr F 2 ) 

105868 v A 60 ; 

506251846536678653 1 
+ 336866181120 ^60 r 



-Tr F 2 ) 7 

.2^6 

■Tr F 2 ) 4 
■Tr F 2 ) 3 



666245292706591 /rTn » Wm nlJ .,1 „ 2 . 4 

Tr F 8 Tr F 12 )( — Tr F 2 4 

2412453120000 v A A 60 7 

575488458230809 , m „ 8wrT , „ Uw 1 m „ 2 . 

(Tr F 8 ) Tr F 14 ) — Tr F 2 

9649812480000 v A A 60 

572769851 f8 F i8 )( l Tr F 2 } 

1766476800 1 A A 60 ; 
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1090453 , m ~,, 

+ Tr F 8 ) Tr F ) 

6133600 v A ; 

20281446070347775183 F 8 )( J_ Ti p2)W 
56144363520000 1 A 60 j 

- 19127927 (Tr F 8 ) 2 (Tr F 12 ) 
2770944000 v ' v 7 



1076029152338172071 , m 1 m 

(Tr F 8 ) 2 (— Tr F ) 

561443635200000 v 7 v 60 ; 

13886644775887 , m 0,0,1 

(Tr F 8 3 — Tr F ) 

36092805120000 v ; v 60 7 

3432221 12 F i4 )( J_ Tr F 2 } 

8131200 v A A 60 ' 

108447101840059177 tm ....1 »,„ 

H Tr F 12 )(— Tr F 2 ) 8 

969830400000 v A 60 ; 

7079594327 F i 2)2( l Tr F 2 )2 

1941730560 v ; v 60 1 

23843919029848581 tm „ Uw 1 m 

Tr F 14 — Tr F ) 

574393600000 v A 60 ; 

157757 . m , 4x2 

+ Tr F 14 ) 2 

1936000 v 1 

1607603265138373 /m 1 m „ 2 .= 

+ Tr F 18 ) — Tr F ) 

344636160000 v A 60 7 

61788727534443 , m 1 m », 4 

(Tr F 20 )( — Tr F 2 4 

54567392000 v A 60 ' 

286680771 , m „ Mw 1 m » x2 

+ Tr F 24 — Tr F 2 2 

6654560 v A 60 ; 

565720659186764144207 1 2 14 

58817904640000 l 60 j 1 j 



B The computer program for the reduction into invariants 

We have implemented the reduction algorithms into a computer program^] in the language of 
FORM This language is particularly suited for these types of problems. Because of all the 

problems with reduction identities when the number of vertices becomes large, we have restricted 
the program to the case of no more than 16 vertices. If the user needs to run the program with 
more vertices, it can be extended by analogy, but many new reduction identities would have to be 
derived. Alternatively one could decide to not reduce a number of contractions with /'s in them 
and leave them for later evaluation. The user should be warned however that some diagrams with 
16 vertices may need quite some computer time and resources for their evaluation. 

The program consists of three parts. The first part reduces all traces of matrices which do not 
belong to the adjoint representation. Much attention is given to a potential contraction of indices. 
The special cases have been written out in one highly nested loop to take the maximum benefit of 
these contractions. This saves much work when we have to use the algorithm of equation ( j30|) for 
the remaining trace. It is quite useful to rewrite each invariant immediately with the notation of 
equation ( |100| ). This removes invariants which have more than one line contracted with the same 

6 This program can be obtained from "http://norma.nikhef.nl/~t68/FORMapplications/Color" 
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Rep. 


Dimension 


i 

J2 

30 


h 


h2 


hi 


(1,0,0,0,0,0,0,0) 


248 


1 


1 


1 


1 


(0,0,0,0,0,0,1,0) 


3875 


25 


-17 


223 


-521 


(2,0,0,0,0,0,0,0) 


27000 


225 


393 


2073 


8961 


(0,1,0,0,0,0,0,0) 


30380 


245 


119 


-1801 


-7945 


(0,0,0,0,0,0,0,1) 


147250 


1425 


-801 


-3921 


90423 


(1,0,0,0,0,0,1,0) 


779247 


8379 


357 


64677 


-207291 


(3,0,0,0,0,0,0,0) 


1763125 


22750 


64330 


653050 


3872050 


(0,0,1,0,0,0,0,0) 


2450240 


29640 


576 


-300624 


-407160 


(1,1,0,0,0,0,0,0) 


4096000 


51200 


59264 


-176896 


-1416448 


(0,0,0,0,0,0,2,0) 


4881384 


65610 


-68202 


1623078 


-5978610 


(0,0,0,0,0,1,0,0) 


6696000 


88200 


-64176 


344544 


2464392 


(1,0,0,0,0,0,0,1) 


26411008 


372736 


12544 


-928256 


20640256 


(2,0,0,0,0,0,1,0) 


70680000 


1083000 


991440 


15398400 


1956600 


(0,1,0,0,0,0,1,0) 


76271625 


1148175 


-64071 


732969 


-67564287 


(4,0,0,0,0,0,0,0) 


79143000 


1404150 


6100842 


97389402 


723747954 


(0,0,0,1,0,0,0,0) 


146325270 


2360085 


-942669 


-20062029 


81822195 


(0,2,0,0,0,0,0,0) 


203205000 


3441375 


3576615 


-53721225 


-402564225 


(2,1,0,0,0,0,0,0) 


281545875 


4843800 


10500696 


50453496 


233862408 


(0,0,0,0,0,0,1,1) 


301694976 


5068800 


-4540800 


36284160 


244435200 


(1,0,1,0,0,0,0,0) 


344452500 


5740875 


3591945 


-51773175 


-133939575 


(1,0,0,0,0,0,2,0) 


820260000 


14773500 


-7295820 


368355300 


-1650963300 


(1,0,0,0,0,1,0,0) 


1094951000 


19426550 


-3552406 


74388314 


342570098 


(2,0,0,0,0,0,0,1) 


2172667860 


40883535 


36197469 


182867949 


3494811285 


(0,1,0,0,0,0,0,1) 


2275896000 


42214200 


-2179296 


-439704336 


1139298552 


(0,0,0,0,0,0,3,0) 


2903770000 


60885500 


-95237740 


3610174100 


-17484769900 


(3,0,0,0,0,0,1,0) 


3929713760 


79228100 


167887580 


2727186380 


10466026340 


(0,0,0,0,0,0,0,2) 


4076399250 


83281275 


-77203203 


-459950403 


17875089309 


(0,0,1,0,0,0,1,0) 


4825673125 


93400125 


-36251565 


-180072525 


-4026565725 


(0,0,0,0,1,0,0,0) 


6899079264 


139094340 


-107301348 


-484327668 


13082745060 



Table 9: Indices for Eg, part one. 



f in a very natural way without any extra pattern matching. For this type of rewriting FORM even 
has a special statement (ToVector). 

The second part eliminates the loops of /'s. Here we do not have to worry about contracted 
indices inside the trace, because FORM looks each time for the smallest loop to make its next 
reduction. Hence this part of the program is much simpler. Because this reduction is much faster 
than the general reduction the first routine calls this reduction routine each time after it has 
removed contracted indices and generated more /'s. Very often such a removal generates a loop 
of /'s and if this is removed it may introduce new contracted indices. The net result can be a 
significant increase in speed. In some cases this is not faster however. Hence there is the option 
not to hunt for /-loops until all other representations have been rewritten. 

The third part of the program contains the reduction identities. Here the program tries to 
eliminate contractions with /'s that are obviously zero and to rewrite the contractions for which it 
can construct meaningful identities. This is a rather peculiar piece of programming. The derivation 
of the reduction identities is by no means a fixed algorithm. The equation (|96| ) can be applied to 
one of the cLr and one of the /'s, but it is not always clear which pair will give good results and 
which will make things worse. In general it seems to be a good strategy to try to increase the 
number of contractions between the invariants (number of common indices). In that case there 
will be more contractions between the /'s and hence more chance of loops that contain only /'s. 
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/l8 


ha 


hi 


^30 


1 


41 


199 


61 


-281 


720023 


-8538743 


107370139 


131601 


20785953 


1677921087 


32641770621 


-130825 


-21485681 


-1669283839 


-32749110565 


-3057657 


1091333799 


69614416281 


-13332825829797 


3122949 


-891843603 


-70053431037 


13392095601009 


158684770 


53898275690 


19220216027590 


4181215457807530 


96664440 


42322995216 


18320504001024 


4178331237939240 


-129144448 


-47849709056 


-18734041431424 


-4186512285307648 


-131658210 


336853672758 


-90053043268518 


53479289344535190 


32750232 


-373796000256 


71666802925296 


-57656847733698648 


-887500544 


601951099904 


-72720785548544 


50171949070921216 


1422281640 


-496637707680 


79406160046320 


-50034441438444840 


368688753 


-283297282191 


-5188200203889 


7356020195237373 


52000832994 


23915855965002 


18712064549246118 


18617881056766529514 


10867212915 


-105084506709 


-9407941453913691 


- 16544807547516665625 


-32456806065 


-12291541050705 


-4580281352771055 


-1096019298618319485 


-14097200472 


-10058442950664 


-14025344636105496 


-17530899917498537112 


-22110000000 


4164533921280 


9773219242968960 


16439802976484140800 


13626625545 


10661778671985 


13958084762321535 


17581170727669309605 


-10809258420 


79696508512740 


-32082728508187860 


-3180131072643324180 


19892530658 


-97084816917686 


17486065830031526 


- 14251422467900989462 


-107677764315 


118967521198869 


-19927893621180789 


12855700625814459585 


-134688633528 


112636179273744 


-29916313535701344 


-3637801450738060968 


-675793438060 


2261894204588980 


-1168369134581957620 


7589921677844796014660 


721988940500 


90443981815180 


94451334234961220 


3296282167094693780 


-1346442094371 


689963459361477 


507354515738178603 


-7162275330089795125911 


518531279235 


147692804386155 


73753893160393605 


17639532662269149015 


104563589460 


-617925299730228 


-535510167643211772 


712889047128020301 1860 



Table 10: Indices for Eg, part two. 



This is not always possible in a direct way, and sometimes we have to just try equation (|9^) in the 
hope that in the next pass the improvement will follow. The selection of the invariant and the / 
that take part in this game has to be done carefully and the code consists of two pieces that make 
a slightly different choice. By running a loop that contains the first choice twice and the second 
choice once, we could reduce nearly everything up to 16 vertices. Of course we had to define the 
object of formula ( |112|) as a separate entity. Similarly we had to define three such objects at level 
16. The program also defines the topologies at level 16 that we did not present in the text. 



C Some examples 

We have run a number of color traces with the program. Here we present the results with some 
timing information. All runs were done on a PentiumPro 200 processor running NextStep. First 
we look at traces of the type 

Tr[T£ ■■■T£T l R 1 ■■■T i g] 

as such traces represent some type of maximal complexity. Here we show the results to n = 7 in 



table 11 



Actually the program can go to n = 8. For this it took about 1520 sec. We do not give the answer 
here. Similarly we can calculate this in the adjoint representation only. This is of course much 
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n 


time 


result 


2 


0.23 s 


N A h{R)(C R -C A /2) 


3 


0.23 s 


N A I 2 (R)(C R - C A ){C R - C A /2) 


4 


0.25 s 


d abcd d abcd + n a I 2 (R)(C r - 3 C 2 R C A + 11/4 C R C' A - 19/24 C\) 


5 


0.95 s 


d R bcd d a A bcd {5C R - 6C A ) + l/3/ 2 (R) d a A bcd d a A bcd 

+N A I 2 {R){C% - 5 C%C A + 35/4 C\C\ - 155/24 C R C\ + 125/72 C 4 ) 


6 


2.59 s 


-8d a R bcdej d A bcde * + Qd abcd d abe} d cdei + I 2 {R)d A be} d cde} {2 C R - 199/60 C A ) 

+d a R bef d c A def (15 C\ - 87/2 C R C A + 179/6 C|) + N A I 2 (R)(C R 

-15/2 C^,Ca + 85/4 C^- 115/4 C|C^ +905/48 C^C^ -1405/288 C^) 


7 


34.9 s 


+112/3 d abcdet d a A bC9 d A }g - 328/9 d a A bc } d a R bC9 d d f 9 
+d R bcdef d a A bcdef (-56 C R + 296/3 CU) 

+< cd 4 6e/ 4 de/ (42 C R - 749/10 C A ) + 67/15 I 2 (R)d a A bcd d a A bef d A ef 
+d abcd d abcd^ C 3 _ 357/2C% C A + 868/3 C fl C A - 2695/18 C\) 
+I 2 {R)d a A bcd d a A bcd (7 C 2 R - 1603/60C fi C A + 497/20 C\) 
+N A I 2 {R){+C% - 21/2 C\C A + 175/4 C%C\ - 280/3 
+5215/48 C\C\ - 19075/288 C R C\ +43357/2592 C\) 



Table 11: Results for traces of the type Tr[T R x ■■■T R n T l R 1 --- T R n }. 



faster because the program selects automatically the smallest loops. These results can be found 
in table For n = 8 the program took 1.5 sec We notice that here the computer time does not 



n 


time 


result 


2 


0.15 s 


±N A C A 


3 


0.20 s 





4 


0.23 s 


jabcdjabcd I at /~i<i 
a A a A ~ 24 iV ^°A 


5 


0.78 s 


\C A d A bcd dT i - ^N A C\ 


6 


0.81 s 


Mbcdjabe} icdef l n2 Mbcdjabcd 13 7\r , r<6 

a A a A a A + jy A a A a A - %^iv A ^ A 


7 


0.89 s 


Srabcdef,abcgrdejg 53 n labcd^abef icdej 
~9 a A a A a A W UA a A a A A 

-fm n aC a 



Table 12: Like the previous table but now in the adjoint representation. 

increase very much with the number of crossing lines (the number of vertices and hence the weight 
is 2n). There is actually more 'compilation time' than 'execution' time. The jump in time going 
from n = 4 to n = 5 represents the use of the reduction algorithms to eliminate /'s. In that case 
the program needs considerably more compilation time. 

The fact that diagrams with only vertices in the adjoint representation are easier to evaluate 
than the diagrams with vertices in the other representations is exactly the opposite of what happens 
with the Cvitanovic algorithms JIJ. For them each / is converted to one or more terms with one 
or more matrices in the fundamental representation. This can lead to an avalanche of terms at the 
intermediate stages, because no advantage is taken from the potentially simpler structures. In the 
case of the traces in the fundamental representation the Cvitanovic algorithms are much faster. 
These algorithms do not worry about symmetrizations and are directly applicable to such traces. 

As an example of high complexity for purely adjoint diagrams we take the one topology of girth 
6 with 14 vertices. It is also called the Coxeter graph. In this the smallest loop has 6 vertices. The 
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result is rather short: 

G 6 (n = 14) 
This took 1.6 sec. 



abode f jabcg rdefg 



raocg,aejg - r MbcdJ 
A d A a A ~ YE A A 



abcd^abef ,cdef 



A d A+^AC 7 A (183) 



D Explicit expressions 

Here we present some expressions for a number of invariants. These are mostly invariants for 
representations that can be used as reference representations. The expressions are given in terms 
of the normalization factor r\ defined in (p7|). 

In all cases the tensors (referred to as d n (R)) are defined as 



d a,...a„ 



Str T^ 1 ... T 1 ^" . 



(184) 



In particular no traces are subtracted and no overall factors are included. For SO(N) we deviate 
from the preferred index normalization of table (|j), since otherwise we would have to deal with a 
few low-iV SO(N) cases separately. 

Results for the fundamental (vector) representation V of SU(N): 



C v 
C A 
d 33 (VV) 

d M (VV) 

d 55 (VV) 

d^VVV) 

dm(VV) 

d 6 3 3 (VVV) 

d bA3 {VVV) 

duiiVVV) 

d3333(VVVV) 



2aN 



1) 



^(iV 2 -l)(iV 2 -4) 

(iV 2 - 1)(N 4 - 6JV 2 + 18) 



6N 2 
a 5 
24N 3 
a 5 



(N 2 - 1)(N 2 - 2) (TV 4 + 24) 



6N 2 

.6 



(TV 2 - l)(iV 2 -4)(iV 2 -6) 
" (iV 2 - l)(N 8 + 6iV 6 - 60iV 4 + 600) 
(Ar2 _ i) 2 (AT 2 _4) 2 

(iV 2 - 1)(N 2 - 4)(iV 4 - 6iV 2 + 18) 



120/V 4 

a 6 
480/V 3 
a 6 



288/V 3 



27N 3 



(N 2 -l){N e 



97V 4 + 817V 2 



189) 



8N 2 



(N 2 - l)(N 2 -4)(N 2 - 12) 



(185) 
(186) 

(187) 
(188) 
(189) 
(190) 
(191) 
(192) 
(193) 
(194) 
(195) 



with a=\r) (see equation (p7|)). The choice rj = l;a= \ corresponds to the most commonly used 
normalization. Then Tr T V T V = \8 ab . In SU(N) the vector representation is always equal to the 
reference representation. 

For the vector representation of SO(N) we have: 



C v 
C A 



a 
2 

a{N 



(N-l) 
-2) 



(196) 
(197) 
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N A = In{N-\) (198) 

d u (VV) = ^N A (N 2 -N + A) (199) 

= ^N A (N A + 2) (200) 
a 6 

dm{VV) = -—N A (N 4 -2N 3 + 33N 2 -32N + 52) (201) 
iy2o 

= ^N A (N 2 A + 16N A + 13) (202) 
a 6 

d444(VVV) = — iV A (2iV 3 - 3N 2 + 33iV - 16) (203) 

(204) 

with a = 7] (see equation (|27|)). In this case a = rj = 2 is the most frequently used convention. Then 
Tr T$T\r = 25 ab . 

Note that for SO(N), N < 6 the vector representation is not the reference representation. The 
formulas for the reference representation for those groups can be read off from the appropriate 
SU(N) or Sp(N) results. The tensors used in the foregoing formulas are normalized so that Ia{Y) = 
I&(y) = 1 in (p5|). As explained above, this differs from the index normalization chosen in table 
(0) for SO (7) and SO (8). For these groups our convention is to make I4 twice as large, and hence 

twice as small. 

For the vector representation of Sp(N) we have: 

Cv = |(iV + l) (205) 
C A = a(N + 2) (206) 
N A = In{N + \) (207) 
a 4 

d u {VV) = —N A {N 2 + N + A) (208) 
a 4 

= — N A (N A + 2) (209) 
a 6 

dmiyV) = -—N A (N 4 + 2N 3 + 33N 2 + 32N + 52) (210) 
1920 



= ^-N A (N 2 A + 16N A + 13) (211) 
a 6 

d4M{VVV) = — N A (2N 3 + 3iV 2 + 33iV + 16) (212) 

(213) 

with a = 7jT\ (see equation (|27])). In this case the vector coincides with the reference representation. 

For all groups for which Ia{A) = we can derive a number of invariants with relatively sim- 
ple methods. This is of particular interest for the exceptional algebras, which have Ia{R) = 
for any representation. We will present the following formulas for tensors defined in the adjoint 
representation, which is not the reference representation (except for E$). The reason for doing 
this is that it allows us to write a single set of relations for all algebras. It is straightforward to 
re-express these results in terms of the reference representation. To do so one needs the relation 
between adjoint tensors and reference tensors, which follows directly from the characters of both 
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representations; the latter can be computed by means of the methods used in Appendix A. The 
adjoint representation is unsuitable for the odd traces of Eq, which are discussed separately below. 
To do the computations, first we notice that 



cf™ 4 = I 22 (A)(5 lll2 5 l3H + S lll3 6 l2H + 5 llH 5 l2ls )/3 



(214) 



\C\, \N A + 2) 



Using this and equation ( |118|) we can determine d§§{AA). Next we 



with I 2 2(A) 

can go even further by using a technique similar to the one used to derive equation ( 118 ): We run 
the program for the product of two traces with 8 vertices. First we run it for a representation R 
in one trace and the adjoint representation in the other. After the run we put R equal to A. This 
gives an expression that includes d$s(AA) and dg^^AAA) and objects that contain combinations 
of d A and d%. We can run the same traces with the program, but starting with both of them in the 
adjoint representation. In that case we obtain an expression that does not contain d%%(AA). This 
gives us the required equation. Now we substitute d A and we need an equation for d% 1 '" 6 which 
is also easy to obtain with the program: 



21 A 



1 



+ ^(d% lP2l3 d% 41516 + 



■)/10 



(215) 



in which we have to take the 10 symmetric combinations over the indices in the last term. For d$ 



we have a similar equation which is given by equation ( 142 ). In total we obtain: 



d M (AA) 
d m (AA) 
dm(AAA) 
deu(AAA) 
d 8 g(AA) 

d 8U (AAA) 

d m4 (AAA) 
d iAAAa {AAAA) 
d i4 44b(AAAA) 



25 C\ 



12(N A + 2 
C\N A ( 797 
' y 288 



(N A + 2)2 
C\N A 



(N A 
175 



27 

,125 125 AT . 



1 

864" 



f 2)2 v 27 
C\N A 



216 



48 (N A + 2) 



3425 
U008 + 



C\N A 
(N A + 2f 

C%N A 125 , 
(N A + 2) 3 1 24 288 A) 



111025 
145152 
625 



N A + 



125 
6804 



N 2 A + 



25 



435456 



C\N A 
{N A + 2Y 

C\N A 
{N A + 2)3 

C\N A 



,5455 3485 ,., 



864 
3125 
324 
6875 



2592"" 2592 

625 Ar 625 

H N A H 

216 1296 



N\) 



3125 »r X 

+ tttttt:^) 



(JVa + 2)3 v 648 ' 1296 
The last two topologies are defined as 

dAteAa(Pl,P2,PZ,Pi) 
d A 4A4b(Pl,P2,P3,P4) 



(216) 

(217) 
(218) 
(219) 
(220) 
(221) 
(222) 
(223) 
(224) 




(225) 
(226) 
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group 


G 2 


F4 




E 7 


Eg 


N A 


14 


52 


78 


133 


248 


C A 


4 T] 


9 77 


24 77 


18 77 


30 77 



Table 13: Values of Na and Ca for the exceptional groups. 



The appropriate values to be substituted for the various groups are given in table 13. 
For Eq we also have to consider the invariants with 5 indices. These we can obtain with the 
reference representation r. For invariants that involve this representation we have: 





52 

" Y V 




N r 


= 27 


(228) 


^44 (rr) 


= 18720 rf 


(229) 


d 55 (rr) 


= 291200 t/ 5 


(230) 


d m (rr) 


= 1023152 r/ 6 


(231) 




= 536640 rf 


(232) 




11211200 , 
3 11 


(233) 


G?644 (rrr) 


20616544 7 
= 5 11 


(234) 


d 554 (rrr) 


= 582400 rj 7 


(235) 


d 88 (rr) 


466346192 „ 
27 ' 


(236) 


of these quantities can be obtained 


in various ways and serve as a check of our pro£ 


rrams. 



E Chiral representations of SO(2m) 

The algebra D m has an index of order m that vanishes for the vector representation, and is non- 
zero only for chiral representations. The latter are characterized by having unequal values for the 
last two Dynkin labels, a m _i 7^ a m . The simplest representations of this type are the spinors, 
S = (0, . . . , 0, 1, 0) and S' = (0, . . . , 0, 1). If m is odd S' is the complex conjugate of S. Note that 
for odd m the extra index (henceforth referred to as the "chiral index") has odd order, unlike all 
other SO(N) indices, whereas for m even there are two distinct indices of even order, namely the 
chiral index and one of the regular indices. 

Since the chiral index vanishes for the vector representation, whereas all indices are non-zero 
for the spinor, it might be argued that the latter is perhaps a more suitable choice for the reference 
representation. However, the vector representation has other advantages, perhaps most importantly 
that the trace-identities for traces of order larger than 2m are simpler. The main drawback of this 
choice is that it requires a separate discussion of chiral traces, which we give here. 

Our conventions for SO (2m) are as follows. The generators in the vector representation are 

M£ = i^fi/2(5?6» - (237) 

This is a complete set of generators for fi < v. The normalization factor \Jnj2 is introduced in 
(p7|). Clearly for SO(N) the most attractive choice is 77 = 2. 



44 



The generators of the (chiral plus anti-chiral) spinor representation are then 

^ = y^j2[^,Y] , (238) 

with {7^,7^} = 2<7 /ii/ . The dimension of this Clifford algebra is 2 m . 
For the chiral 7-matrix (often referred to as 75) we choose 

7c =(i)V---7 2m (239) 

The phase is chosen so that 7 C is hermitean (assuming the other 7-matrices are chosen hermitean 
as well). Hermiticity fixes the overall factor up to an m-dependent sign, for which we have made a 
conventional choice. 

The chiral spinor generators are 

£f = ^(l± 7c ) (240) 
The symmetrized trace of m such generators is 

StrE^ 1 . . . T,^ Vm = ^Strl!^ 1 . . . ^ mUm ± ^StrE^ 1 . . . £^ m 7c (241) 

The first term yields (for m even) an ordinary d-tensor, the second one yields the d tensor. By 
general arguments, 

StrE^ 1 . . . Yl^ Vn = I n (S±)d^ lv ^-^ y ^ + lower order 

±I n {S±)^ v ^-^ u ^ , (242) 

where [/^i^i] denotes an adjoint index. By convention 

I n (S±) = ±1 (243) 

then 

tf[p-\v\},...,\p im v m \ _ IgtrJ]' 111 ' 1 . . . Y} lmVm ^ c = -(rj/2) m / 2 e ,J ' lr "' l ' m , (244) 

with the definition g 1 ----' 2 " 1 = ]_. Note that the interpretation requires pairs of indices to be identified 
with an adjoint index, and that the e tensor is indeed fully symmetric under permutation of pairs 
of indices. 

This gives us an explicit expression for the extra tensor, and computing invariants that involve 
this tensor is now straightforward. One can either do that by computing all other tensors also 
in the spinor representation. Then computing the invariants amounts to simple 7-algebra. Since 
expressions for the spinor characters in terms of vector traces are explicitly known (see (^5|) and 
(pig)), all tensors in the spinor representation can be re-expanded in terms of reference tensors. 
Alternatively it is also very easy to work directly with the reference tensors, which, using (|237| ) 
can be expressed completely in terms of kronecker <5's with vector indices. These indices are to be 
taken in pairs and identified with adjoint indices, and can then be pairwise contracted with the e 
tensor. 

For the normalization of the tensor d of SO '(2m) we get 

d mm (CC) = ^2" m (2m)! 
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The factor 2 m compensates the double counting of index pairs [iv as compared to adjoint indices. 
The argument "CC" indicates that two chiral tensors d are used. Furthermore we have 

dmm(CV) = d rnrn (CS) = . 

Here "V" is the fundamental representation, as before, and "S" refers to the non-chiral tensor 
computed in the spinor representation, 

S 



2 



This tensor is related to the one in the reference representation V. 

To order 12, the maximal order we used for the other SO(N) tensors, the only other quantity 
of interest is ^444, which has chiral tensor contributions only for SO (8). Here is a complete set of 
results for that group. 



d AA (VV) 


= 70r? 4 


(245) 


d AA (SS) 


8 


(246) 


d AA {CC) 


_ 315 4 
— 8 


(247) 


dm(VV) 


_ 581 6 
~~ 8 '/ 


(248) 


d AAA {VVV) 


= 70t] 6 


(249) 


d AAA {SSS) 


_ 665^6 
— 32 '/ 


(250) 


d AAA (SSC) = d AAA (VVC) 


= 


(251) 


d AAA (SCC) 


_ 525 Jo 
~ 32 '/ 


(252) 


d AAA (VCC) = 


_ -105 6 
8 


(253) 


d AAA {CCC) 


= 


(254) 



The fundamental quantities here are the ones involving d A (V), de(V) and d A {C). The tensor d A (S) 
is related to d A (V): 



jabcd jabcd , ° jabcd 

a S — ~o% + o fl 2,2 



1 jabcd , 3 

2 dy + 2 

The results involving d A (S) are given here because one may use them to check the triality relations 



d AA (VV) = d AA (SS) + 2d AAA (SC) + d AA {CC) 

and 

d AAA {VVV) = d AAA {SSS) + 3d AAA (SSC) + 3d AAA {SCC) + d AAA {CCC) 

which are due to the fact that the tensors d A (V) and the combinations d A (S) ± d A (C) are given by 
symmetrized traces of the triality-related representations (8„), (8 S ) and (8 C ). 
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